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1.9. Literature citations 


References are normally cited in the text by placing the name(s) and 
the year, without any comma between them, in parentheses. If there 
are two authors for one citation, both names should be given, sep- 
arated by an ampersand (&). If there are more than two authors, only 
the first name should be given, followed by “et al.”. Commas should 
only be used to separate two or more years linked with the same author 
(author group). If two or more citations are made in one set of 
parentheses, they should be separated by a semi-colon. If more than 
one citation fora particular author (author group) is made for the same 
year, “a”, “b”, “c”, etc. should be added to the year. If citations are 
made within the normal running text, only the year(s) should be 
placed in parentheses. The following examples will illustrate the 
required style: 
— (Copernicus 1986) 
— (Copernicus & Galilei 1988) 
— (Hubble et al. 1985; Newton et al. 1987; Ptolemaus & Copernicus 
1989a,b, 1992) 
— Recently Galilei et al. (1992, 1993) showed that... 
Authors’ initials are permitted only in exceptional cases; for example, 
to distinguish between two authors with the same surname. Each 
literature citation made in the text should have a corresponding entry 
in the References at the end of the paper (see Sect. 1.13 below). For 
frequently cited papers an abbreviated form of citation is recom- 
mended, e.g. Paper I, Paper II (if appropriate) or by the initial letters 
of the authors’ surnames. The abbreviation should be defined at the 
first place of mention in the text. 
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1.13. References 


The reference list should contain all the citations occurring in the text, 
ordered alphabetically by surname (with initials following). If there 
are several references to the same first author, they should be entered 
according to the following scheme: 

(a) One author: chronologically. 

(b) Author + one co-author: alphabetically by co-author, then chro- 
nologicaily. 

(c) Author + two or more co-authors: chronologically. 

The number of authors listed in a single reference should be limited to 
five. If a particular work was written by more than five authors, only 
the first three should be listed, followed by “et al.”. This rule also 
applies to the number of editors in the case of edited works. 

In the new styling for the references, introduced in 1990, the use of 
italics and boldface has been abolished, and the punctuation has been 
simplified. Simplified abbreviations are now used for the commonest 
astronomy journals. A list of these is given in Sect. 1.13.1 below. 
Standard abbreviations should be used for other journals, following 
the recommendations given by the IAU (see Sect. 1.13.2 in the 
unabridged version of the /nstructions). Asample reference list show- 
ing the required style is given in Sect. 1.13.3. 


Journals: References to journals should normally only contain the 
authors’ names, the year. the name of the journal (abbreviated), the 
volume number and the number of the first page. Neither the title of 
the paper itself nor the last page number should be given. 


Monographs: An entry for a monograph should include the title of the 
book, the name of the publisher, the place of publication and, if 
appropriate, the relevant page number(s). 


Edited works: Entries for handbooks, series, conference proceedings, 
etc. should include the title of the relevant chapter, the name(s) of the 
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editor(s), the title of the complete work, the name of the publisher, the 
place of publication and the first page number. 

Items which are not yet published may be included in the list of 
references and denoted by “in preparation”, “submitted”, or “in 
press”. As faras possible authors should update these references when 
making their corrections in the page proofs. This also applies to items 
characterized as “preprint” or “private communication” if these have 
turned into publications at a later date. 


1.13.1. Simplified abbreviations of frequently used journals 


ARA&A Annual Review of Astronomy and Astrophysics 

Afz Astrofizica 

AJ Astronomical Journal (the) 

AZh Astronomiceskij Zhurnal 

A&A Astronomy and Astrophysics (Letters indicated by the 
page number) 

Astronomy and Astrophysics Supplement Series 

Astronomy and Astrophysics Review (the) 

Astrophysical Journal (the) (Letters indicated by the 
page number) 

Astrophysical Journal Supplement Series (the) 

Astrophysics and Space Science 

Bulletin of the American Astronomical Society 

Journal of Astrophysics and Astronomy 

Monthly Notices of the Royal Astronomical Society 

Nature 

Publications of the Astronomical Society of Japan 

Publications of the Astronomical Society of the Pacific 

PASPC PASP, Conference Proceedings 

QJRAS Quarterly Journal of the Royal Astronomical Society 

Sci Science 

SvA Soviet Astronomy 


A&AS 
A&AR 
ApJ 


ApJS 
Ap&SS 
BAAS 
JA&A 
MNRAS 
Nat 
PASJ 
PASP 
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Sample reference list 
Attention should be paid to the order of the items in each reference 
and to the punctuation used. 


Bohr N., Einstein A., Fermi E., 1992, MNRAS 301, 257 (BEF) 

Curie M., Curie P., 1991, A&A 248, 612 

de Gaulle C., 1996, Solar Astronomy. Oxford Univ. Press, Oxford 

Heisenberg W., Yang C.N., 1993, Chin. Astron. 537, 36 (Paper III) 

Laurel S., Hardy O., 1994, Active Galactic Nuclei. In. Churchill W., 
Roosevelt F.D., Stalin J. (eds.) The Evolution and Distribution 
of Galaxies. Wiley, New York, p. 210 

Shakespeare W., Goethe J.W. von, Moliére J.-B., et al., 1995, Super- 
nova 1994D. In: Duck D., Mouse M. (eds.) Proc. [AU Symp. 408, 
The Physics of Supernovae. Reidel, Dordrecht (in press) 


4. Letters to the Editor 
4.1. General remarks 


Every attempt will be made to ensure that Letters are published within 
6-10 weeks after receipt of the camera-ready manuscript by the 
publisher. Initially 3 copies of a draft version of the Letter should be 
submitted to the Letter Editor. He will then contact the authors 
regarding provisional acceptance and also send full details on the 
procedure to be followed for the submission of the final camera- 
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ready version. A few essential points may be mentioned here. Letters 
are limited to four printed pages. The publisher requires a camera- 
ready manuscript of good quality. Diskettes should not normally be 
submitted. Manuscripts in which the text or the artwork is of inad- 
equate quality will not be accepted for publication. 


4.2. Title page 


Since stylistic corrections to the title, authors’ names and addresses 
are often necessary, and since the received/accepted line usually still 
has to be inserted, the title block will always be typeset separately and 
later combined with the camera-ready text. Thus the text for the title 
block should always be typed or printed on a separate sheet. 


4.3. Camera-ready pages 


Asarule, a vertical space of 8—9 cm should be left free on the first of the 
camera-ready pages so that the typeset title block can be inserted. 
However, the exact amount of space needed depends to some extent 
on the length of the title and the number of addresses. It is recom- 
mended therefore that authors who use one of the Springer-TpX 
macro packages (P-AA vers. 3.0 or L-AA vers. 3.0) for writing their 
text include the title block in the printout of the first page. In this way, 
the correct page layout will be achieved, even though the title block 
will be replaced later by the typeset version. 

It is recommended that authors use one of the two Springer-TpX 
macro packages for generating their text since the required page 
layout and print quality will be achieved without much difficulty. 
Details on how to obtain these packages are given on the first page of 
these /nstructions. 

Alternatively, authors may wish to work with a word processing 
program (such as Word or WordPerfect) in combination with a page- 
layout program (such as Pagemaker or Ventura). In this case, the 
output must have proportional letter spacing and the exact col- 
umn/page format of A & A. The typographical style of A & A should 
be followed as closely as possible. 

Line drawings should be reduced to the required size and pasted 
into the text at the appropriate positions. This also applies to tables if 
these have been written separately. The empty frames generated by 
the Springer-T—X macro packages to mark the sizes and positions of 
figures should be covered or removed when the figures are pasted in. 
Line drawings, which are submitted separately and have not been 
reduced to the final size cause extra work and increase the risk of 
delay. However, halftone and coloured illustrations are an exception to 
this rule: they should not be pasted into the text because they will have 
to be screened, reproduced separately and inserted later when the film 
pages are made up. 


Using one of the Springer-Tp, X macro packages, authors will have 
a choice of producing a text output in either PostScript (PS) Times 
fonts or in Computer Modern (CM) fonts. However, the PS output is 
slightly more compact than the CM output and this fact should be 
taken into account if the length of the text is critical with respect to the 
page limit prescribed for Letters. 


5. Permission to reprint 


Requests for permission to reprint figures or tables which have already 
been published in A.& A. should be addressed to the copyright 
holders, not to the publisher. The address to write to is as follows: 
Professor J. Lequeux, Astronomy and Astrophysics Editorial Office, 
Observatoire de Meudon, F-92195 Meudon Principal Cedex, France. 
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Abstract. We present a multifrequency statistical analysis of 
spiral galaxies belonging to the Virgo cluster using photomet- 
rical data available in the literature or taken by us. The analysis 
shows that in spiral galaxies the global star formation activity 
is tightly coupled with the atomic gas content, while it is al- 
most independent of the molecular gas reservoir (as determined 
from the CO lines). In spiral discs the surface brightness of stars 
emitting in the visible and in the infrared H band is correlated 
with the H2 column density. These observational evidences are 
interpreted as a result of using the CO emission of spiral galax- 
ies to infer their global molecular gas content. This technique, 
in fact, is valid only in solar-neighbourhood environments, with 
abundant molecular gas and star formation, while it leads to an 
underestimate of the H2 content in those regions of spiral discs 
which have abundant cold molecular gas and low mass-star for- 
mation activity. In high-mass star formation regions, instead, 
the intense far-UV radiation field produced by young O-B stars 
partially photodissociates the diffuse molecular gas, weakening 
the expected relation between star formation and molecular gas 
content. 

In the Virgo spiral galaxies the molecular gas content as 
traced by CO is not influenced by the cluster environment. Even 
galaxies showing strong signs of interaction with the intergalac- 
tic medium, as determined from their perturbed H I morphology, 
have a normal molecular gas content. 


Key words: galaxies: spiral — galaxies: ISM — galaxies: clus- 
ters: Virgo cluster — galaxies: intergalactic medium — galaxies: 
interactions 





1. Introduction 


The comparative analysis of a multifrequency data set is an 
appropriate tool for understanding the role played by the in- 
terstellar gas on the activity and on the evolution of galaxies. 
Observations of the H 1 line at 21 cm and of the various CO lines 
at millimetric wavelengths give information on the total content 
and on the distribution of the atomic and molecular gas, nec- 
essary to form stars (Young & Scoville 1991). Utraviolet and 


optical data give clues on the activity of star formation in galax- 
ies (Kennicutt 1989). Radio continuum emissivity (synchrotron 
radiation of relativistic electrons accelerated in supernovae rem- 
mants) (Condon 1992; Gavazzi & Jaffe 1986) and FIR luminos- 
ity (emission of the interstellar dust heated by stars) are indirect 
indicators of the star formation activity in spiral discs. 

The analysis of the relationships between these different 
emissions provides information on the mechanisms that regulate 
the transformation of gas into stars in spiral galaxies. Further- 
more such studies can help us to understand how these mech- 
anisms are affected by the environmental effects, i.e. by ram 
pressure sweeping and tidal interactions. 

Multifrequency analysis of different samples of galaxies 
have led to different and sometimes controversial results. In 
spiral galaxies the activity of star formation is strongly corre- 
lated with the HI content (Scodeggio & Gavazzi 1992; Kenni- 
cutt 1989) but, contrary to what would be expected, it is only 
weakly correlated with the molecular gas reservoir (Kennicutt 
1989; Buat 1992). Cluster galaxies are in general H I-deficient 
(Chamaraux et al. 1980; Haynes et al. 1984), indicating that 
the ram pressure exerted by the intergalactic medium can effi- 
ciently remove the neutral hydrogen. Conversely the molecular 
gas distribution is not influenced by an interaction with the dense 
intracluster gas, as expected in a ram pressure scenario (Kenney 
& Young 1989; Casoli et al. 1991; Young & Scoville 1991). This 
statement, however, has been recently criticized by Rengarajan 
& Iyengar (1992), who found that the molecular gas content of 
Virgo cluster galaxies is on average smaller than that of field 
objects. Furthermore some spiral galaxies in the Fornax cluster 
have been found with a normal HI content but deficient in Hz 
(Casoli et al. 1994). Preliminary results, however, of CO emis- 
sion in an unbiased survey (Casoli et al., in preparation) suggest 
that the FCRAO sample of Young et al. (1989), used as refer- 
ence by Rengarajan & Iyengar (1992) and Casoli et al. (1994), 
is CO-rich compared to average galaxies, then indicating that 
Virgo and Fornax cluster galaxies have a normal molecular gas 
content. The activity of star formation, instead, does not seem 
influenced by the cluster environment. Cluster galaxies have on 
average star formation rates similar (Gavazzi et al. 1991; Donas 
et al. 1990; Kennicutt et al. 1984) or marginally lower (Kenni- 
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Table 1. The available photometric data for the Virgo spirals (number of galaxies) 





Spiral H: CO F100 F60 H 


J-K B-V U-B_ UV24 UV20 UVIS5S 


Ha 





151 145 75 151 151 86 54 
Detect 129 62 118 113 


108 80 


75 31 33 37 
65 28 22 





cutt 1983a; Kodaira et al. 1990; Bicay & Giovanelli 1987) than 
that of field objects. Furthermore, some cluster galaxies have 
been observed with deficient H1 content but with a high activ- 
ity of star formation, which is probably induced by the recent 
interaction of these galaxies with the intracluster gas (Bothun 
& Dressler 1986; Gavazzi 1989; Gavazzi et al. 1992; Boselli et 
al. 1994a; Caldwell et al. 1993). 

The aim of this work is to analyze the relations between the 
gas content and the activity of star formation in spiral galaxies 
and to understand if these mechanisms are influenced by the en- 
vironment where galaxies evolve. We analyze a sample of spiral 
galaxies belonging to the Virgo cluster (Sect. 2) with multifre- 
quency data available in the literature or taken by us (Boselli et 
al. 1994b). In Sect. 3.1 we analyze the available star formation 
indices. In Sect. 3.2 we discuss the relations between star for- 
mation activity and gas content, both atomic and molecular, in 
spiral galaxies. In Sect. 4 we interpret the new results. The ef- 
fects of environment on the gas content of galaxies are analyzed 
in Sect. 5, and general conclusions are reported in Sect. 6. 


2. The sample 


The sample analyzed in this paper is taken from the CGCG cat- 
alogue (Zwicky et al. 1961-1968) limited to the photographic 
magnitude mp, < 14.5 in the region containing the Virgo clus- 
ter (12" < a < 13" and0° < 6 < 18°) and with morphological 
classification later or equal to Sa. The sample consists of 151 ob- 
jects. Due to the relatively bright photographic magnitude limit 
adopted, no irregular and dwarf galaxies are included in the sam- 
ple. The morphological classification is taken from Binggeli et 
al. (1985) when available, or from de Vaucouleurs et al. (1991; 
RC3). We use the same membership assignmentents as adopted 
by Binggeli et al. (1985). We assume a distance of 11 Mpc for 
Virgo members, 23 Mpc for possible members, M and W clouds, 
14 Mpc for the southern extension and the distance determined 
from the redshift (with Hyp = 100 kms~! Mpc!) for back- 
ground objects in order to take in account the 3-dimensional 
distribution of galaxies inside the cluster. In particular the dis- 
tance for possible members is assumed to be 23 Mpc since these 
galaxies have on average a higher recessional velocity with re- 
spect to that of the cluster members and are mainly located in 
the western and southern part of the cluster, then closer to the 
higher velocity M and W clouds and to the southern extension 
(Binggeli et al. 1993). 

All galaxies have redshift available, and the majority of them 
have been measured in the 21cm HI line and in the far in- 
frared (Table 1). Half of the sample has data available in the CO 
line, optical colours, near infrared magnitudes and UV fluxes at 


2421 A. Ha and UV fluxes at 2000 and 1560 A are available for 
a few galaxies (Table 1). 

H1 data are from Warmels (1986), Cayatte et al. (1990), 
Bottinelli et al. (1990), Schneider et al. (1990) and references 
therein. CO fluxes are taken from Kenney & Young (1988), 
Stark et al. (1986) and Boselli et al. (1994b), and references 
therein. Far infrared data are from Lonsdale et al. (1985), Soifer 
et al. (1989), Haynes et al. (1990) and Condon & Broeder- 
ick (1991). Near infrared and optical data are from Boselli & 
Gavazzi (1994) and from Longo & de Vaucouleurs (1983) and 
references therein. UV fluxes are from Smith & Cornett (1982; 
at 2421 A), Donas et al. (1987; at 2000 A) and Kodaira et al. 
(1990; at 1560 A). Ha data are taken from Kennicutt & Kent 
(1983) and Romanishin (1990). 

Hydrogen molecular masses have been determined adopting 
the data reduction procedure described in Boselli et al. (1994b). 
We used a standard conversion factor of (Strong et al. 1988): 


N(H2) = 2.3 107° Ico(K kms~') in molem~? 

H>2 masses relative to each single beam have been determined 
using the relation: 

M(H2) = 1.06 1077 (dist/1 pe)? 2.3Ic¢0(@/1')?> in Mo 

where N(H2) (in mol cm~*) is the molecular hydrogen column 
density, Ico (in K km s~') is the intensity of the 12CO(1-0) line, 
and @ is the HPBW of the used telescope. Galaxies with large 
angular dimensions compared to the beam size of the telescope 
were mapped in contiguous beams (Kenney & Young 1988; 
Stark et al. 1986) and the total mass M(H2)r is simply given 
by the sum of the 1/(H2); determined in each single position i, 
then by the relation: 
M(H2)r = YM (H2); for 1 = ie n. 

Only the positions with a detection are used in order not to 
add unwanted noise. Galaxies with smaller angular dimension 
have been observed only in the central position (Boselli et al. 
1994b). Since galaxies are undersampled, lower limits to the 
global H2 masses are obtained. However, as shown by Boselli 
et al. (1994b), since the molecular gas is generally centrally 
peaked (Kenney & Young 1988), total Hz masses should not 
differ significantly from the single-position estimate. 

Total HI and H2 masses have been estimated adopting the 
average distance of the subgroup of the Virgo cluster to which 
each galaxy belongs, according to Binggeli et al. (1985). Opti- 
cal magnitudes are corrected for Galactic extinction (Burstein & 
Heiles 1982), internal extinction (Boselli & Gavazzi 1994) and 
K-dimming (de Vaucouleurs et al. 1976, RC2). Near infrared 
magnitudes are corrected for internal extinction as described in 
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where a and b are the optical UGC diameters (Nilson 1973), 
and adopting D = 1.27 irrespective of the morphological clas- 
sification. This correction is then a usual conversion of the UV 
magnitudes to face-on orientation. We preferred this correction 
technique to the one proposed by Donas & Deharveng (1984), 
based on the hydrogen content of each object, because their tech- 
nique would introduce a systematic and artificial difference be- 
tween galaxies with a normal HI content and H I-deficient ones. 
The H1 gas located outside the optical disc, in fact, depleted in 
HI deficient galaxies, would only marginally contribute to the 
extinction of the UV fluxes. 

In the following analysis we will often compare the Virgo 
galaxies with the larger sample of relatively isolated galaxies 
0.4 0.8 12 1.6 in the Coma supercluster and with objects in other clusters de- 

Log Halpha E.W. (A) scribed in Gavazzi et al. (1991; hereafter GBK) to analyze the 
effects of the environment on galaxy evolution. Unfortunately 
these larger databases do not contain CO and UV measurements. 


-2 


UV(2421 A) surface brightness 
mag arcmin 








| | i J 








Table 2. Linear regression parameters between the different star for- 
mation indices. The following regressions are determined using only 
detected objects. When a « is marked, the regression have been deter- 
mined using the sample of spiral galaxies described in GBK 
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mag arcmin 


DUV24, = —1.35(+0.38) log Ha + 16.63(+0.60) R? = 0.45 obj = 17 
DUV24, = 2.97(+0.48)(U — B), + 15.27(+0.77) R? = 0.53 obj = 36 
DUV24, = 4.01(+0.69)(B — V). + 13.12(+0.83) R? = 0.41 obj = 50 
DUV20, = —1.61(+0.47) log Ha + 15.92(+0.44) R? = 0.49 obj = 14 
DUV15, = —2.28(+0.42) log Ha + 17.85(+0.33) R? = 0.83 obj =8 
FIR/H = 0.75(+0.07) log Ha + 0.32(+0.32) R? = 0.48 obj = 122* 
FIR60/a? = 0.63(+0.10) log Ha — 1.14(+0.46) R? = 0.24 obj = 131* 
FIR100/a? = 0.56(+0.09) log Ha — 0.60(+0.39) R? = 0.26 obj = 120* 


UV(2421 A) surface brightness 

















3. Discussion 
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3.1. Star formation indices 


Many indicators of the star formation activity (SFA) of galaxies 
have been proposed in the literature (see the review by Kennicutt 
1990). The Ha equivalent width (E.W.), which is probably the 
most widely used, is an index of the recent (< 10’ years), high 
mass (M > 10M,) star formation activity (SFA) of galaxies 
(Kennicutt 1983b; 1990). The paucity of data at this wavelength, 
0.2 0.4 0.6 0.8 however, induced many authors to search for other star forma- 
c (B-V)c tion indices. (B — V), and (U — B), corrected colours, which 
Fig. la—c. Relations between the different optical and ultraviolet ete well conelated ye te Ho EW. (Rennsoutt 1983b; GBK), 
star formation indices: UV surface brightnesses at 2421 A (in mag mies indicators of the integrated star formation history of galax- 
arcmin~”) versus a log Ha E.W. (in A), b(U — B), ande(B—V), _1€S (even if they also depend on the present day star formation 
rate averaged over the last 108 years) (Larson & Tinsley 1978; 
Lequeux 1988). UV surface brightnesses are good star forma- 
Boselli & Gavazzi (1994). UV observed magnitudes are cor- tion indices of the young stellar component of intermediate mass 
rected for internal extinction using the relation: (O’Connell 1990; Buat et al. 1987). In Fig. 1 we plot the hybrid 
UV corrected surface brightness DUV,(2421A), defined as: 


UV(2421 A) surface brightness 
mag arcmin 











Am, = —2.5D log(b/a) SUV, = mag UV, + 5 log a mag arcmin™” (1) 
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where a is the UGC optical diameter, versus log Ha E.W., 
(U — B), and (B — V),. The regression parameters between 
the different star formation indices are listed in Table 2. UV 
surface brightnesses are good star formation indices, being well 
correlated with Ha and with optical colours. The corrected UV 
surface brightness (Eq. 1) correlates better with the Ha E.W. at 
shorter wavelengths (see Table 2): this indicates that the UV sur- 
face brightness at 1560 A is a good indicator of the recent high 
massive star formation, while at longer wavelengths (2421 A) 
the UV surface brightness is also sensitive to the less recent 
star formation history of galaxies. Together the Ha E.W., UV 
surface brightness and optical colours can be used as indicators 
of the global SFA of stars of different mass and age. Corrected 
optical and, more severely, UV surface brightnesses may how- 
ever still be affected by extinction. Other indirect star forma- 
tion indices, which have the advantage of being not affected by 
extinction, have been proposed in the literature, such as nor- 
malized FIR fluxes (at 60 and 100 microns; Persson & Helou 
1987; Devereux & Young 1990; Young et al. 1989) and radio 
continuum fluxes (Kennicutt 1990; Buat 1992). The FIR flux 
(at 60 zm) normalized to the near infrared H band flux (de- 
fined as in GBK) and the FIR hybrid surface brightnesses at 60 
and 100 ym are plotted in Fig.2 versus log Ha E.W. FIR/H is 
reasonably well correlated with the Ha E.W. (GBK) and then 
it can be used as a star formation index. Only a trend with a 
large scatter is present between the FIR hybrid surface bright- 
nesses (both at 60 and 100 zm) and the Ha E.W. (see Table 2), 
indicating that FIR surface brightnesses are not good star forma- 
tion indices, this result being confirmed by the analysis of some 
nearby galaxies resolved in the FIR (Rand et al. 1992 and refer- 
ences therein). For this reason the FIR surface brightnesses will 
not be used as star formation indicators. Indeed many authors 
(Persson and Helou 1987; Boulanger & Pérault 1988; GBK; 
Sauvage & Thuan, 1992; Rand et al. 1992) pointed out that at 
FIR wavelengths the global emission of galaxies is not only due 
to dust heated by young stars, but also to heating by the general 
interstellar field. 

Synchrotron emission, finally, may be contaminated by ther- 
mal emission and by nuclear activity. 

In conclusion the Ha E.W. is an indicator of the high massive 
SFA (O-B stars) of spiral galaxies during the last 10’ years. The 
UV surface brightness is an indicator of the formation of stars 
of intermediate mass over time scales of about 10° years; at 
short wavelengths (1560 A) the UV luminosity traces the SFA 
of stars of higher mass over shorter time scales than at longer 
wavelengths (2421 A). Optical (U — B), and (B — V), colours 
are indicators of the integrated star formation history of galaxies 
(< 10° years). Since these star formation indices are affected 
by extinction, we will also use throughout this work the FIR/H 
ratio, with the caveat that FIR luminosities are related not only 
with the emission of dust heated in the H II regions by O-B stars, 
but also by the diffuse light from all stars. 
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Fig. 2a—c. Relations between a the FIR/H ratio, the logarithm of the hy- 
brid FIR surface brightnesses at b 60 and c at 100 um (in Jy arcmin~*) 
and the log Ha E.W. (in A). The linear regression parameters are listed 
in Table 2 


3.2. Star formation and gas surface densities 


Neutral hydrogen is the primary gaseous component neces- 
sary to form new stars through its condensation into molecu- 
lar clouds. Early studies (Schmidt 1959) proposed that the star 
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Table 3. Coefficients of the linear regression analysis between the the star formation indices and the log of the H1 and the H2 surface densities 
and the H1 gas mass (normalized to the dynamical mass). Only detected galaxies are included 





log Ha = 0.58(+0.09)SH1 — 2.59(+0.22) 
(U — B). = —0.25(+0.05)5H1 + 1.71(+0.21) 
(B — V)- = —0.15(+0.03)SH1 + 1.54(+0.15) 
FIR/H = 0.57(+0.08)5H1 — 2.65(+0.33) 


log Ha = —0.03(+0.15)SH2 + 1.40(+0.34) 
(U — B), = —0.01(+0.08)"H2 + 1.12(+0.23) 
(B — V), = 0.10(+0.04)SH2 + 0.01(+0.14) 
FIR/H = 0.24(40.14)SH2 — 0.45(+0.44) 


log Ha = 0.49(+0.08) log(Mu1/Mwx) + 2.37(+0.21) 


R? = 0.57 
R? =0.31 
R? =0.22 
R? = 0.45 
R? = 0.00 
R? = 0.00 
R? =0.10 
R? = 0.06 


R? = 0.62 


obj=37 
obj=69 
obj=94 
obj=71 
obj=27 
obj=43 
obj=55 
obj=48 
obj=27 
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Fig. 3a—d. Relations between the star formation indices a Ha E.W. (in A), b (U — B),, ¢ (B — V), and d FIR/H and the logarithm of the HI 
surface density (in Me kpe~”). The linear regression parameters are listed in Table 3 


formation rate (SFR) in galaxies is related with the local gas 
surface density of atomic hydrogen by a simple power low of 
the type: 


SFR =@2° . 


Kennicutt (1989), analyzing the relation between star formation 
and gas surface densities in resolved galaxies, concluded that 
this relation is correct only when high values of gas surface 


densities are present, while for gas surface densities lower than 
a certain threshold, star formation does not take place, or in 
other words that the relation between the gas content and the 
star formation rate is non-linear. Adopting the available star 
formation indices discussed in the previous section, we analyze 
the relationship between the SFA and the atomic and molecular 
gas content in Virgo spiral galaxies. Since not all galaxies are 
resolved in HI and CO, we use hybrid surface densities of the 
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Fig. 4a-d. Relations between the star formation indices a Ha E.W. (in A), b (U — B)., ¢ (B — V), and d FIR/H and the logarithm of the H2 
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surface density (in Me kpc~“). The linear regression parameters are listed in Table 3 


gas content defined as: 


Degas = log( Meas /d”) 


Mo kpe~? 


where d is the linear diameter determined from the optical UGC 
major diameter; galaxies distances are determined as described 
in Sect. 2. Since optical diameters are in general larger than CO 
diameters (Kenney & Young 1988) and smaller than the HI 
ones (Cayatte et al. 1990), molecular and atomic gas densities 
determined with this relation are respectively under-and over- 
estimated. As the morphological distribution of the atomic and 
molecular gas component are quite different from each other, 
the total hybrid H 1+H) surface density does not have a physical 
meaning and will not be used in this work. Scodeggio & Gavazzi 
(1992), comparing optical and HI radii of a sample of field and 
cluster resolved galaxies, concluded that the error introduced 
by using hybrid H1 gas densities is small and it can in practice 
be neglected. 

The various star formation indices are plotted in Figs. 3 and 
4 versus the hybrid surface densities SH 1 and SH) (in logarith- 
mic scale). The coefficients of the different linear regressions 
are listed in Table 3. The SFA is tightly correlated with the 


Table 4. Results of the linear regression analysis between the stellar sur- 
face brightnesses and the gas surface densities (in logarithmic scales). 
Only detected galaxies are included 





LA, = 0.26(+0.17)5H1 + 9.81(+0.82) R? = 0.03 obj = 82 
ov = —0.11(+0.11)0H1 + 14.57(+0.54) R? = 0.01 obj = 94 
UB, = —0.29(£0.08)H 1 + 16.35(+0.48) R? = 0.09 obj = 128 
SU, = —0.53(40.12)5H1 + 17.92(+0.52) R? = 0.23 obj = 69 
DUV24. = —1.20(+0.19)OH1 + 23.02(+0.74) R? = 0.44 obj = 55 
DUV20, = —1.48(+£0.30)SH1 + 23.74(40.61) R? = 0.48 obj = 28 
SUV 15. = —1.59(+0.31)OH1 + 24.98(40.67) R? = 0.58 obj = 21 


LH, = —1.00(+0.19)H2 + 17.40(+0.65) R? = 0.35 obj = 53 
xv = —0.78(+0.11)SH2 + 18.56(+0.38) R? = 0.49 obj = 55 
XB, = —0.68(+0.10)0H2 + 18.54(+0.37) R? = 0.43 obj = 62 
“SU, = —0.73(+0.18)XH2 + 18.86(+0.51) R? = 0.30 obj = 43 
DUV24, = —0.23(+0.36)OH2 + 16.87(+1.01) R? =0.01 obj = 31 
DUV20, = —0.49(+0.56)ZH2 + 16.87(+0.81) R? = 0.05 obj = 16 
DUV 15. = 0.79(+0.64)OH2 + 10.81(40.99) R? =0.13 obj = 12 





HI! surface density: galaxies have higher star formation when 
HI is higher, i.e. when there is more HI available (Scodeg- 
gio & Gavazzi 1992; Buat 1992; Kennicutt 1989; Guiderdoni & 
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Fig. 5a and b. Relations between the SFA and the gas mass normalized 
to the dynamical mass: a log Ha E.W. (in A) versus log MH1/Myx, b 
log Ha E.W. (in A) versus log MH2/My. The linear regression pa- 
rameters are listed in Table 3 


Rocca-Volmerange 1985). Notice that because of the adopted 
definition for the gas surface density, this relation indicates that 
the SFA is correlated not only with the HI gas located over 
the optical disc of each galaxy, as shown by Kennicutt 1989, 
but also with the one located outside the main star forming re- 
gions, in principle physically disconnected to the SFR. Part of 
the spread in SH 1 (Fig. 3) is due to the existence of many HI- 
deficient galaxies in the sample (Scodeggio & Gavazzi 1992). 
The Hp» surface density, instead, does not correlate with any 
star formation index (only very marginally with FIR/H). Buat 
(1992) and Kennicutt (1989), analyzing the relations between 
the gas surface density and the hybrid UV2000 and Ha surface 
brightness on smaller samples of spiral galaxies, reached similar 
conclusions. 

The Ha E.W. is plotted in Fig.5 versus the atomic and 
molecular gas masses normalized to the dynamical mass. 
Galaxy total dynamical masses (M,o,) are determined using the 
virial theorem as described by Gavazzi (1993), Miot = AVE? 
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Fig. 6a and b. Relations between the surface brightness of the stellar 
component emitting a in the UV at 1560A and b in the V band (in 
mag arcmin~”) and the logarithm of the hybrid HI surface density (in 
Mo kpc~’). The linear regression parameters are listed in Table 4 


(for galaxies with inclination > 30° and good H 1 quality spec- 
tra), where AVyjc is the deprojected 21 cm line width and r 
the optical radius. Notice that M,. is an hybrid entry since it 
is determined with H1 widths and optical radii (see Gavazzi 
1993). The linear regression parameters are listed in Table 3. 
Figure 5 indicates that the activity of star formation traced by 
Ha E.W. is strongly correlated with MH /M,, while it is not 
with My, /Meor. 

The comparison of the H1 and Hp» surface densities with 
the surface brightness of the different stellar component may 
provide informations on the way gas condenses to form stars. 
The hybrid stellar surface brightnesses in UV at 1560 A and in 
V (corrected for internal extinction) are plotted in Figs. 6 and 
7 respectively against the HI and H2 surface densities. The re- 
sults of the linear regression analysis are listed in Table 4. As 
expected, young and massive stars are more frequent in galaxies 
with high H I surface densities; the relation between the star for- 
mation indicators and 5H 1 is stronger in Ha and UV (Scodeggio 
& Gavazzi 1992; Kennicutt 1989; Buat 1992; Buat et al. 1989) 
than in the optical (B, V) and near infrared (H). The surface 
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Fig. 7a and b. Relations between the surface brightness of the stellar 
component emitting in a the UV at 1560 A and b in the V band (in 
mag arcmin~*) and the logarithm of the hybrid H2 surface density (in 
Mo kpe~?). The linear regression parameters are listed in Table 4 


brightness of the young stellar component (Ha, UV) instead, is 
not correlated with the molecular gas surface density, (Kennicutt 
1989; Buat 1992; Buat et al. 1989), while the surface brightness 
of the stellar component emitting in the optical (V, B, U) and 
NIR (#) bands correlates with the SH) (notice that in the rela- 
tion ©V — SHzp the scatter is quite low). 


4. Interpretation 


A surprising result of the present analysis is the lack of corre- 
lation between the SFA and the molecular gas content of spiral 
galaxies as traced from CO. Kennicutt (1989) discussed the pos- 
sible observational uncertainties which could affect this result. 
He concluded that it may be due to (i) the effect of extinction 
on the different star formation indices, (ii) large variations of 
the CO-H>2 conversion factor among different galaxies, (iii) a 
tighter relation of star formation with the gas volume density 
than with the gas surface density, or (iv) the fact that H2 is cou- 
pled with a low mass star formation mode, rather than with the 
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Fig. 8. Relation between the logarithm of the hybrid H2 surface density 
(in Me kpe~?) and the index of metallicity J — K 


massive stellar population traced by Ha. (i) Even if the extinc- 
tion is not negligible, the first hypothesis can be in principle 
ruled out since extinction would introduce the same scatter in 
the relation between star formation and © H 1 (Fig. 3). Further- 
more only the H!I surface density (and not that of H2) is related 
with the star formation index FIR/H, which is not affected by 
extinction. (ii) The CO-H>2 conversion factor may depend on 
the galaxy metallicity (Z) (Kenney 1987; Buat 1992). Since 
in our sample magellanic irregulars and dwarf galaxies are ex- 
cluded (see Sect. 2), the analyzed objects are expected to have 
similar metallicities (Pagel & Edmunds 1981). Unfortunately 
we have direct measurements of metallicity for only 4 galaxies 
in our sample (Shields et al. 1991). Therefore, as suggested by 
Bothun et al. (1984), we use the near infrared colour J — K 
as an indirect index of metallicity. Using the empirical relation 
reported in Thuan (1985) between J — K and the metallicity Z, 
we estimate that the observed galaxies have metallicity close to 
the galactic one (0.25Z) < Z < 2Z ). Furthermore no correla- 
tion is observed between /H)> and J — K (Fig. 8) (the error in 
the J — K colour index is 0.1 mag). All these evidences suggest 
that the lack of relation between the molecular gas content and 
the SFA is not induced by the use of the same (mean) CO-H2 
conversion factor for all galaxies. (iii) Figure 5 indicates that, 
even if the SFA was strictly governed by the normalized gas 
mass rather than by the gas surface density, in spiral discs the 
activity of star formation is correlated with the neutral hydrogen 
content and not to the molecular gas reservoir. (iv) The relation 
between the H2 surface density and the stellar component emit- 
ting in the V, B and marginally in the H band may confirm the 
hypothesis of Kennicutt (1989) that the molecular gas is cou- 
pled with a low mass star formation mode, rather than with the 
massive stellar population traced by Ha. 

Our analysis shows that the global star formation activity of 
spiral galaxies is correlated with the HI gas content and not to 
the molecular gas content as traced by CO. The first fact is not 
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too surprising as it is HI which condenses to form molecular 
clouds and then stars (Seiden 1983). However the second fact is 
more surprising just because in the Milky Way stars appear to 
form in molecular clouds (Shu et al. 1987). In this respect, one 
should first discuss to what extent CO traces molecular hydro- 
gen. While the usual CO line intensity-H2 column density rela- 
tion is valid for the conditions of the solar neighbourhood where 
it has actually been determined (Strong et al. 1988), there are 
circumstances in which it is grossly invalid. Allen & Lequeux 
(1993) have found giant molecular clouds with a CO emission 
10 times lower than that of the Galactic clouds of the same size 
in the inner regions of M31. They interpreted this as due to a lack 
of far-UV radiation from young stars, such that while neither 
CO nor H2 are photodissociated (the observed amount of atomic 
gas in these regions is very small), CO is at a very low kinetic 
and excitation temperature and emits very little. Note that mas- 
sive stars do not form from these clouds: this is a clear example 
of a lack of relation between molecular gas and SFA. Another 
extreme example (in the opposite direction) is that of molecular 
clouds in the actively star-forming regions of the SMC (Rubio et 
al. 1993). Here the far-UV radiation field is very intense due to 
the high density of young stars and the small absorption by dust, 
so that CO and H)2 survive photodissociation only in the dense 
clumps inside relatively small (10-30 pc) clouds. Any diffuse 
gas outside this clouds and even in any low-density medium in- 
side them is photodissociated (Shaya & Federman 1987). As a 
result the CO line intensity-H2 surface density conversion fac- 
tor is different from Galactic, although in the clouds the CO 
line emission is only 2 times weaker than in the correspond- 
ing Galactic clouds, because the smaller surface filling factor of 
the CO-emitting clumps is partly compensated by their higher 
CO emission due to a higher temperature (Rubio et al. 1993; 
Wolfire et al. 1993). Note that in these regions there is active 
star formation with little molecular gas, but abundant H 1, once 
again illustrating the lack of relation of SFA with molecular 
gas or with CO emission but its strong correlation with atomic 
gas. We do not understand why the molecular clouds in M31 do 
not form massive stars (although they may form low-mass stars 
like the cold Taurus clouds in our Galaxy: these stars cannot 
be observed in M31). Another unanswered question is why the 
SFA is so active in the inner regions of the SMC where there is 
so little molecular gas (Rubio et al. 1991) (although it remains 
clear that molecular gas, while not abundant, is a necessary step 
in star formation). 


It is likely that galaxies contain at the same time three 
regimes: a solar neighbourhood-like one, corresponding to 
abundant molecular gas with “normal” CO emission and star 
formation, a SMC-like one with little molecular gas and CO 
emission but very active massive star formation, and a regime 
with abundant cold molecular gas, low CO emission and weak 
massive-star formation. An interesting fact is that while a large 
fraction of molecular gas can be in the latter regime as a reser- 
voir for (perhaps) future massive-star formation, it contains little 
atomic gas; atomic gas is abundant in the other regimes repre- 
senting 1/2 of the total gas in the solar neighbourhood, and 
almost 100% in active star forming regions where there is little 


molecular gas (as observed in M51 and M83: Allen et al. 1986; 
Vogel et al. 1988; Tilanus & Allen 1989, 1991; Rand & Kulkarni 
1990; Rand et al. 1992). CO traces mainly the moderately ac- 
tive star-forming regions. The interplay between these regimes 
and the reason why some regions containing much H) (but with 
little CO emission) do not form stars is not understood, but the 
Ico to N(H2) conversion might be strongly affected in a way 
which would blur out a global relation between CO and SFA. 


5. Environment 


A study of the environmental effects on the HI content in clus- 
ter galaxies has been made by Haynes et al. (1984) using a 
well-defined reference sample of isolated objects. Unfortunately 
no unbiased comparison samples of isolated galaxies with CO 
and Ha data are yet available forbidding similar studies on 
the molecular gas content and on the SFA of spirals. The CO 
FCRAO survey (Young et al. 1989), often used as a reference 
sample, is in fact biased towards bright infrared galaxies and in- 
teracting systems. The sample of field galaxies of Gavazzi et al. 
(1991), instead, is not complete in Ha. Since no reference sam- 
ples of isolated galaxies with both CO and Ha data are available, 
we analyse if inside the Virgo cluster, as for the H 1 component, 
the molecular gas reservoir as traced by CO changes with the 
projected distance from the cluster centre. The H1 deficiency 
parameter (defined as in Haynes & Giovanelli 1984) can also 
be used as an index of the efficiency of the interaction between 
the galaxies and the intracluster medium: in galaxies with an 
HI! deficiency > 0.3 the H1 gas has been efficiently removed 
during the interaction. The H! deficiency parameter is plotted 
versus the angular distance to the X-ray centre of the cluster in 
Fig. 9. These two quantities are well correlated, indicating that 
Virgo galaxies closer than 5° from the X-ray centre are gener- 
ally H1 deficient, and that there are no objects with a normal H! 
content at a projected angular distance < 2° from the cluster 
centre (Chamaraux et al. 1980; Warmels 1986; Cayatte et al. 
1990; 1994 and references therein). 

The Hp surface brightness is plotted in Fig. 10 (a) versus 
the angular distance from the cluster centre and (b) versus the 
H 1 deficiency parameter. The lack of any relation indicates that 
the molecular gas content (as traced by CO) of Virgo galax- 
ies is not influenced by the environment. Many Virgo objects 
have been resolved in HI: Cayatte et al. (1994) divided spi- 
ral galaxies into 4 different groups according to their HI gas 
morphology. Group I includes objects with an H1 distribution 
similar to the one found in field galaxies, with a ratio of the 
HI to optical diameters Dy1/Dopn > 1.3 (Dy obtained at 
an HI isophote of 107° at cm~?). Group II includes galaxies 
with 0.75 < Dy1/Dopt < 1.3, group III includes galaxies with 
Dyi/Dopn < 9.75. Galaxies belonging to group IV (anemic; 
0.6 < Dy1/Dop < 0.85) have H1 distributed on a low surface 
density ring with a strong central hole. Cayatte et al. (1994) in- 
terpreted the group membership as an indicator of the duration 
and of the strength of the interaction between the galaxy and 
the intracluster medium. Group I members did not have enough 
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Fig. 9. Relation between the H I-deficiency parameter and the angu- 
lar distance from the X-ray center of the Virgo cluster (in degrees). 
Filled circles are H I detected galaxies, while empty triangles represent 
undetected objects 


time to interact with the intergalactic medium, while in group 
IV the interaction was so efficient that HI has been completely 
removed. 

Group IV members have normal H2 hybrid surface bright- 
nesses: DH) = 5.74; 5.84; 5.90 Mokpe~? for NGC 4450, 
4548, 4579 respectively compared with the average value for 
the Virgo cluster sample (only detected galaxies) of SH2 = 
5.99 + 0.46Mo kpe~?. This indicates that the molecular gas 
content is not influenced even in the most perturbed objects of 
the sample (Boselli et al. 1994a). 

Stark et al. (1986), and Kenney & Young (1989), reached 
the same conclusion analyzing a smaller sample of Virgo spiral 
galaxies. A similar result was found by Casoli et al. (1991) 
on the Coma cluster, where even stronger interactions between 
galaxies and the intragalactic medium are expected. 

These observational evidences are currently interpreted as- 
suming that ram pressure stripping takes place in clusters of 
galaxies (Gunn & Gott 1972). The ram pressure pV? exerted on 
galactic discs moving at V ~ 1000kms~' through the inter- 
galactic medium is sufficient to overcome the restoring gravita- 
tional force which keeps the H | anchored to the discs of galaxies, 
while it is not sufficient to remove the molecular gas since Hz 
(as traced by CO) is distributed deeper inside the galactic po- 
tential well (Young & Scoville 1991). However the abundant 
cold molecular gas which may exist in the outer regions of spi- 
ral galaxies and normally does not emit in CO (Lequeux et al 
1993) might be stripped out together with H1. 


6. Conclusions 


In this paper we analyze the relations between the activity of star 
formation and the atomic and molecular gas content of spiral 
galaxies and the effect of the environment on the gas content 
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Fig. 10a and b. Relation between the logarithm of the hybrid H2 surface 
density (in Mo kpc”) and athe angular distance from the cluster center 
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of discs using a multifrequency statistical analysis based on 
a sample of 151 Virgo cluster spiral galaxies with available 
photometrical data at different wavelengths. 

We confirm that in spiral galaxies the UV surface brightness 
is an indicator of the formation of stars of intermediate mass over 
time scales of about 10° years. The correlations among the Ha 
E.W., (U — B), and (B—V), optical colours and the UV surface 
brightnesses at various wavelengths indicate that at 1560 A UV 
luminosities trace the SFA of stars of higher mass and over 
shorter time scales than at longer wavelengths (2421 A). The 
FIR surface brightnesses, instead, are only weakly correlated 
with the SFA of spiral galaxies. 

This analysis shows that the SFA of spiral discs is corre- 
lated with their HI gas reservoir, while it is almost independent 
of their Hz gas content as traced by CO. The stellar compo- 
nent emitting in the optical and near infrared bands, instead, is 
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correlated with the surface density of the molecular gas reser- 
voir. A possible interpretation of these observational evidences 
may be problems in converting the CO line emission to the 
global molecular hydrogen content of spiral discs: the usually 
adopted relation is valid only in specifics environments. Spiral 
galaxies, in fact, may contain simultaneously three regimes: a 
solar-neighbourhood like one, with abundant molecular gas and 
star formation (where the CO emission is a good indicator of 
the H2 gas content), a SMC-like one with little molecular gas 
and CO emission but very active massive-star formation (where 
both H2 and CO have been photodissociated by the far-UV radi- 
ation field), and a regime with abundant cold molecular gas (not 
well traced by the low CO emission) with weak massive-star 
formation activity. 

Contrary to what happens to the atomic gas reservoir, the 
molecular gas content of the Virgo spiral galaxies (as traced by 
CO) is not influenced by the hostile cluster environment. Even 
galaxies with stronger signs of interaction, as determined by 
their H! morphological distributions, have normal H) contents. 
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Abstract. In an effort to explain the ionization of gas in elliptical 
galaxies, we have computed the Lyman continuum output from 
an ageing starburst, using isochrone synthesis spectral evolu- 
tion techniques. We have shown that post-AGB stars provide 
sufficient ionizing radiation to account for the observed Ha lu- 
minosity and equivalent width in very early type galaxies. Pho- 
toionization models using the radiation field from these stars are 
shown to easily reproduce the excitation level of the ionized gas. 
Thus, we argue that the ionization observed in normal elliptical 
galaxies is provided by their old stellar populations. 


Key words: galaxies: elliptical — galaxies: ISM — galaxies: stel- 
lar content 





1. Introduction 


Elliptical galaxies, once thought to be devoid of gas, are now 
known to possess a quite diversified interstellar medium (ISM). 
Although the detection rate of H I remains deceptive (3 out of 64 
ellipticals in the Bregman et al. 1992 sample), evidence of cold 
gas is provided by thermal dust emission in the infrared. As- 
suming a canonical dust-to-gas ratio of 10~7, the total amount 
of cold gas inferred is 10°—10° Mo. The mass of gas in the hot 
phase as determined by X-rays is in the range 10’-10'! Mowith 
a detection rate around 70% (Bregman et al. 1992 and refer- 
ences therein). The existence of a warm ionized gas component 
is revealed by emission line measurements. The spectroscopic 
survey of Phillips et al. (1986; hereafter PJDSB) reported a 55— 
60% detection rate of [N11] in both ellipticals and SO’s. The 
amount of warm ionized gas inferred is typically in the range 
10°—10* assuming a density of 10° cm~*. A new analysis of 
the PJDSB data by Buson et al. (1993) showed the existence 
of a direct proportionality between the emission line luminosity 
and the integrated blue luminosity (Lem « aaa ). In their 
work on imaging 15 early-type galaxies (mostly ellipticals) with 
spatially extended emission lines, Buson et al. showed that the 
warm ionized component is probably closely related to the cold 
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component of the ISM, and may represent only about 1% of its 
mass. 


The question of the origin of this gas is not yet settled. 
It could be either external to the galaxy, resulting, for exam- 
ple, from a cooling flow or due to accretion from a companion 
galaxy, or internal to the galaxy through stellar mass loss. Breg- 
man et al. (1992) suggest that the hot and cold forms of the 
ISM may have separate origins. Kim (1989) and Bertola et al. 
(1992) argue that the ionized gas present in elliptical galaxies 
is in the majority of cases of external origin. Determination of 
the chemical composition of this gas may give some clues to its 
origin. 

The chemical composition, though, can be safely deter- 
mined only if ionization is understood. There are principally 
four plausible causes for the ionization. One is linked to the 
observed X-ray gas. While recombination in the hot phase as 
inferred from the cooling flow rate fails to account for the Ha 
luminosity by two orders of magnitude in early type galaxies 
(Kim 1989; Trinchieri & di Serego Alighieri 1991), heat transfer 
from the hot gas to cooler gas (Sparks et al. 1989) or photoion- 
ization by the X-rays radiated by the cooling hot gas may be 
sufficient to account for the extended emission lines (Voit & 
Donahue 1990; Donahue & Voit 1991; but see Kim 1989). A 
second possibility is that the ionization is related to some kind 
of nuclear activity (Kim 1989; Shields 1992). As shown by 
PJDSB, the extended emission line regions in elliptical galax- 
ies bear many similarities with LINERs (low ionization nuclear 
emission regions; Heckman 1980) which are considered a pow- 
ered down version of the Seyfert phenomenon. Most models 
of LINERs assume that ionization is provided by an energetic 
radiation originating from the nucleus, usually in the form of a 
power law although other energy distributions have been pro- 
posed as well (Ferland & Netzer 1983; Halpern & Steiner 1983; 
Stasiriska 1984; Péquignot 1984; Binette 1985). A third possi- 
bility is that of shock waves (Heckman et al. 1989). The fourth 
possibility is that the ionizing radiation is arising from some 
component of the stellar population. Independently of the ori- 
gin of the nuclear gas (whether internal or external), we con- 
sider the direct proportionality between the emission lines and 
the blue luminosity as very suggestive of a stellar origin for 





L. Binette et al.: Photoionization in elliptical galaxies by old stars 





T Perey ' T ll T T A | T 


aN 
On 





. © LL | 


—_ Old Post—AGB Stars 
__- Nuclei of young PN 


Young Massive Stars 


4 


Pewee. eo 





-1,,-1 
log Q, (s M ) 
+S o > 
rrr eee eT ee ee eee. Se ©) ee 


nes 
— 


aes 
oO 








| pol 


ul. 1 


1 1 posi 1 


bo aiiiil L 


ae es ae 


Fig. 1. Temporal evolution of the logarithm of Qu 
which is the total number of ionizing photons 
above 13.6eV produced by different subsets of 
stars in an ageing starburst of total initial stellar 
mass normalized to 1 Mo. The dotted line repre- 
sents young massive stars, the dashed line corre- 
sponds to stars which have spent less than 5000 
years in the post-AGB phase and approximately 
represent the nuclei of young planetary nebulae. 
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the ionizing photons. Buson et al. (1993) suggest that the com- 
ponent yielding the ionizing radiation may be the same as the 
one producing the UV upturn observed at 2000 A (Greggio & 
Renzini 1990). Various possibilities have been proposed for the 
ionizing radiation: the ultraviolet emission from hot horizontal- 
branch stars (Minkowski & Osterbrock 1959; Osterbrock 1960), 
from OB stars (Kim 1989; Shields 1992), from evolved massive 
stars (Terlevich & Meinick 1985), and, from post-AGB stars (di 
Serego Alighieri et al. 1990). 

In the present paper, we explore in more details the hypoth- 
esis of photoionization by post-AGB stars. Taking advantage 
of the spectral synthesis techniques described by Bruzual & 
Charlot (1993), we follow the time evolution of a starburst, and 
derive for different epochs not only the number of ionizing pho- 
tons but also their spectral energy distribution. This allows to 
estimate the expected Ha luminosity and the equivalent width, 
assuming complete absorption of the ionizing photons by the 
gas. Finally, using photoionization models, we compute the line 
intensities emitted by the different ions of this gas for various 
mean distances of the emitting regions with respect to the ion- 
izing stars. Comparisons are made with the sample of elliptical 
galaxies studied by Buson et al. (1993) and Demoulin-Ulrich et 
al. (1984). 


2. The Lyman continuum output from an ageing starburst 


We consider a population of stars formed in an instantaneous 
burst of star formation, and distributed following the Salpeter 


The continuous line represents all the older post- 
AGB stars whose ionizing photons merge into the 
diffuse ionizing field 


10°" 


(1955) initial mass function in the mass range between 0.1 and 
125 Mo. The time evolution of the spectral energy distribution 
of this population is computed using the Bruzual & Charlot 
(1993) isochrone synthesis spectral evolution algorithm. The 
library of the stellar tracks used to follow the evolution of in- 
termediate and low mass stars (Mzams < 7 Mo) through the 
post-AGB phase was assembled using the Schénberner (1981, 
1983) and Wood & Faulkner (1986) tracks and adopting the 
initial-final mass relation of Weidemann (1987), as can be seen 
in detail in Magris & Bruzual (1993). The stellar spectra library 
used in this work is basically the same as in Bruzual & Charlot 
(1993). The spectral energy distribution of stars with effective 
temperatures between 25 000 K and 50000 K is taken from the 
Kurucz (1979) model atmospheres. For hotter stars, we have 
adopted a blackbody energy distribution (we shall come back 
to this point below). 


Figure | shows the time variation of Quy, the total num- 
ber of photons above 13.6 eV, emitted by this stellar population 
(normalised to a total initial stellar mass of 1 Mo). For illus- 
trative purposes, we have divided the stellar radiation field into 
three components. The young massive stars are represented by 
the dotted line. The post-AGB stars which have left the AGB 
since less than 5000 yr are represented by the dashed line. In a 
first approximation, these represent the nuclei of young plane- 
tary nebulae. In this case, the ionizing radiation is absorbed by 
the dense planetary nebula shells (as is shown by evolutionary 
photoionization models of planetary nebulae, Stasifiska 1989). 
Finally, stars with post-AGB ages above 5000 yr are represented 
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by the continuous line. Most of their ionizing radiation is leak- 
ing out from the old planetary nebulae shells and is available to 
ionize the general ISM. 

From Fig. | we see that, after about 6 10’ years, the produc- 
tion of ionizing radiation by massive stars is overtaken by that 
coming from post-AGB stars. The latter declines very gently 
during the following 10'° years, remaining at a level of about 
10~> times the maximum attained by the massive stars. Figure 
| also shows that, in an integrated population of stars, the old 
post-AGB stars dominate over the nuclei of young planetary 
nebulae in the production of the Lyman continuum flux. Nu- 
clei of planetary nebulae become completely negligible after 
3 10° years because then stars which have not yet left the AGB 
all have small core masses and require longer than 5000 years 
to reach a temperature necessary for a significant emission of 
Lyman continuum photons (Schénberner 1983). 

Figure 2 shows the spectral energy distribution of the in- 
tegrated post-AGB components, at various epochs. Figure 2a 
gives the energy distribution of the photons emitted by the nu- 
clei of young planetary nebulae at an epoch corresponding to 10? 
years after the burst. For comparison purposes, we have super- 
imposed blackbody energy distributions with effective tempera- 
tures of 7.5 104 K, 10° K and 1.5 10° K. Figure 2b represents the 
photon energy distribution of the old post-AGB stars, at differ- 
ent epochs (10°, 8 10°, 10'° and 1.3 10!° yr after the burst). The 
ionizing radiation from these stars is not trapped by a planetary 
nebula gas shell but rather escapes into the general ISM. We see 
that, as could be expected, the radiation field from the integrated 
post-AGB populations is quite hard, and is not blackbody-like. 
Note however that the amount of photons emitted above 54.4 eV 
is uncertain. Indeed, especially in this domain, the true spectral 
energy distributions for individual post-AGB stars may deviate 
significantly from the blackbody approximation which we have 
adopted in the computations. For example, it highly depends on 
the helium atmospheric abundance, which, for these stars, may 
range from 0.1 to almost one. Also, one should not forget that 
the recent non-LTE model atmospheres, which include spheri- 
cal extension and the effects of stellar winds, still have to face 
the problem of the Zanstra discrepancy in planetary nebulae, as 
discussed in Gabler et al. (1991). However, the uncertainty on 
the radiation field above 54.4 eV will only affect the predicted 
intensities of the He Il recombination lines (cf. Sect. 4) and has 
a negligible effect on the lines arising from ions with lower 
ionization potentials. 


3. The Ha luminosity and equivalent width due to 
ionization by post-AGB stars 


From the values of Qy presented above, one can deduce the 
Ha luminosity resulting from photoionization by old post-AGB 
stars in an ageing starburst, assuming complete reprocessing of 
the ionizing photons. Figure | implies that the ionizing ability 
of an old starburst never declines to zero, within the age of the 
Universe. Thus, a low-level emission line activity in elliptical 
galaxies is a natural consequence of the evolution of the initial 
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Fig. 2. a Integrated energy distribution in »F, (arbitrary units) at a 
burst age of 10° years for stars which have left the AGB track for less 
than 5000 years. For comparison, different blackbody curves have been 
superimposed. b Integrated energy distributions in vF, at burst ages 
of 10’, 8 10°, 10'° and 1.3 10"° years for stars on their post-AGB track 
since more than 5000 years 


starburst, provided cold gas is present. We also note that the 
very slow decrease of Qy due to post-AGB stars, as compared 
with the fast drop during the OB star phase between 10°* and 
10’* years, means that even if ellipticals form initially through 
merging, the Qy corresponding to the current epoch and due to 
evolved stars is simply proportional to the total stellar mass ac- 
cumulated by the merged galaxy. If merging went on for epochs 
exceeding ~ 10° years, however, a precise determination of 
Qu would require taking into account the spread in age of the 
different merging events. Finally, it is obvious that in galactic 
nuclei where star formation occurs today, the contribution of OB 
stars to Qy (and therefore to the luminosity in Ha) would easily 
dominate that of post-AGB stars whenever the stellar mass of 
the starburst exceeds 10~> the mass of all stars contained in the 
nuclear regions. This is consistent with observing much larger 
equivalent widths of Ha in late type galaxies where important 
starburst activity is clearly going on in many of these (Kennicutt 
1992 and references therein; Keel 1983a). 
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Table 1. Observational characteristics of the sample and model predictions 





Object z Ms 2rem/Ae EWnn 


Tem 


(kpc) 


[NuJ/Ha _ Stellar 


Qu 
(pred.) 


L(Ha) 
(pred.) 


L(Ha) 


mass (pred.)/(obs.) 





NGC 1395 —5 
NGC 4636 —5 
NGC 5846 —5 
IC 1459 —5 
NGC 6868 —5 
NGC 1453 —5 
NGC 7097 —5 
NGC 2974 —5 
NGC 3962 —5 
NGC 7302 —3 
NGC 1947 —3 
NGC 484 —3 
ESO 234-G21 —2 
ESO 118-G34 —2 
NGC 745 —1.4 


—21.8 1.1 
—21.3 1.9 
21.9 2.7 
—21.4 1.2 
—21.9 2.0 
—21.9 2.6 
—20.7 1.3 
—21.2 2.0 
—21.2 5.0 
—20.4 0.97 
—-19.9 1.1 
—21.7 4.2 
—21.1 4.9 
—17.8 0.93 
—21.0 2.8 


0.15 

0.20 

0.24 

0.24 

0.25 

0.30 - 
0.32 3.66 
0.44 - 
0.80 _ 
0.31 - 
0.35 2.28 
0.56 1.71 
0.49 15.1 
1.64 14.9 
0.25 2.17 


1.74° 
1.35° 
1.80° 
1.81° 
2.29" 
1.81° 
0.81" 
2.53” 
2.20° 
0.66° 
0.85" 
1.69" 
0.84" 
0.17" 
1.16° 


6.9 10° 
5.6 10° 
1.210 
7.6 10° 
1.210 
1.510 
4.910” 
1.010” 
1.610” 
3.510” 
2.710” 
2.010% 
1.010% 
1.110° 
5.110” 


5.110"! 
4.1 10°! 
8.8 10°! 
5.6 10°! 
9.2 10°! 
1.1 10 
3.6 10°! 
7.6 10°! 
1.2 10° 
2.6 10°! 
2.0 10°! 
1.5107 
7.410"! 
8.210% 
3.8 10°! 


7.010'° 
5.710" 
1.210"! 
7.710" 
1.310"! 
1.510"! 
5.010'° 
1.010"! 
1.710!! 
3.610'° 
2710" 
2.0 10'! 
1.010"! 
1.110" 
5.210'° 


3.46 
1.01 
1.11 
0.43 
1.32 
0.35 
0.59 
0.96 
1.94 
0.53 
0.78 
0.27 
0.04 
0.03 
0.19 





* From Tables VII or II of PJDSB. 
» From Demoulin-Ulrich, Butcher & Boksenberg (1984). 
© From fit of [N n]/Ha vs. Mg correlation in PJDSB (see Sect. 3). 


We now compare the Ha emission expected from the post- 
AGB phase with that measured in the sample of (mostly) el- 
liptical galaxies observed by Buson et al. (1993). Table | gives 
the relevant information for these galaxies: column | lists the 
galaxy names, column 2 their galaxy types (T’) according to the 


Third Reference Catalogue of Bright Galaxies (de Vaucouleurs 
et al. 1991: RC3), column 3 their total absolute magnitude Mj" 
derived from the total blue magnitude By in RC3 assuming 
Ho = 55kms~' Mpc! (adopted throughout this paper), col- 
umn 4 the radius rem (in kpc) of the area in which line emis- 
sion is detected by Buson et al. (1993), column 5, the value 
Of 2rem/Ae, where A, is the effective diameter of the galaxy 
(in kpc) derived from the RC3 catalogue. The ordering in the 
table is in increasing galaxy type (T'), and for objects with the 
same T,, in order of increasing value of 27em/A,. For objects 
in common with the spectroscopic survey of PJDSB, we list 
in column 6 the [N I] equivalent width determined by PJDSB. 
In column 7 we give the [N I1]/Ha ratio measured by PJDSB 
or (in the case of two ellipticals) by Demoulin-Ulrich et al. 
(1984). For the objects for which Ha was too weak to be de- 
tected, we used a simple fit of the correlation found by PJDSB 
between the [N 1!]/Ha ratio and the total absolute magnitude 
Mg": (Nu/Ha = —14.2 — 0.73( Mp" + 5 log 55/ Ho) valid for 
the range —19 > My" > —23. Column 8 is the estimated 
stellar mass inside the region of radius r., using an integrated 
stellar mass-to-light ratio M/ Lx of 8 (consistent with our star- 
burst model at an age of 1.3 10'° yr) and estimating Lg within 
Tem USing the integrated B-band standard aperture-magnitude 
curves of RC3 (Table 11 in RC3). Column 9 lists the number 
of ionizing photons corresponding to the stellar mass within 
the radius rem, using a typical specific photon luminosity of 
7.3 10” quanta/s/Mo(cf. Fig. 1). Column 10 gives the predicted 
Ha luminosity assuming recombination case B and complete re- 


processing of the Lyman continuum photons. It is also supposed 
that partial obscuration by dust lanes does not affect drastically 
the observed Ha flux (we shall come back to this below). Col- 
umn || gives the total Ha luminosity observed inside the radius 
Tem by Buson et al. (after correcting for the blend with [N I] us- 
ing the [N 11]/Ha ratios of column 8). The last column (12) gives 
the ratio of Ha expected over the value observed. 


We find remarkable that for ellipticals most ratios in Column 
12 are of order unity. A significant shortage of ionizing photons 
in the models is only apparent for later types (T. > —4) where we 
might conjecture that some level of stellar (or nuclear) activity 
is currently taking place. It is noteworthy that ESO 118-G34, 
for instance, has a significantly lower [N t1]/Ha (= 0.17) ratio 
which suggests the presence of an HII region (cf. Veilleux & 
Osterbrock 1987). In our example, the adopted age of 1.3 10!° 
years for the starburst is simply consistent with the ages inferred 
by Magris & Bruzual (1993) in their fit of the energy distribution 
of the elliptical galaxies NGC 4472 and NGC 3379. An earlier 
age would in general simply yield more ionizing photons as 
indicated by Fig. 1. 


We have shown that post-AGB stars produce a sufficient 
number of ionizing photons within the region of detected line 
emission. In our model, we interpret the observed steeper radial 
decline in the Ha+[N I] surface brightness compared to the stel- 
lar continuum surface brightness (Buson et al. 1993) as being 
the result of a drastic drop in gas content near radius Tem. In 
this case, no significant reprocessing of Lyman photon would 
be expected to take place at larger radii even though further 
ionizing photons are generated there. A small fraction of these 
photons would be absorbed by the gas rich inner regions, an ef- 
fect not considered here but which would nevertheless enhance 
somewhat the integrated Qy within rem. 
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We have assumed so far that the cold gas intercepts all of the 
ionizing photons produced within rem. If this were not the case, 
an additional source of ionizing photons would be required to 
produce the observed Ha emission. One further complication 
could result if some of the emission line clouds were not visible 
to us owing to large opacity of the cold ISM along the line of 
sight. If this effect were important, it would require from the 
model a production of ionizing photons larger by possibly up to 
a factor two. 

The fact that the emission line gas distribution is far from 
filling the whole region of radius rem in some objects of Table 
| cannot be taken into account in the simple scheme presented 
here. We can, however, compute independently the theoretical 
equivalent width of the Ha emission implied by our starburst 
model! and compare it with values measured in different ob- 
jects and in different positions. Assuming again that all of the 
ionizing photons are absorbed, we find that the Ha equivalent 
width (against the integrated continuum from all the stars con- 
sidered in the burst model) increases from 0.6 A at a burst age of 
8 10° years, to 1.3 A at 10'° years and to 1.7 A at 1.3 10! years. 
A value of 1.7A is certainly consistent with the observations 
of PJDSB, which yield typical values of 1.54 A for the [N 11] 
6583 A equivalent widths which translates into Ha equiva- 
lent widths of 1-3 A (assuming for instance their mean value of 
[N 11]/Ha of 1.4). This independent check on the soundness of 
our model confirms that post-AGB stars are serious candidates 
for generating the amount of ionizing photons implied by the 
warm gas detected in elliptical galaxies. 


4. Photoionization models 


In order to predict the emission line ratios implied by our sce- 
nario, we have computed photoionization models using the pho- 
ton energy distribution corresponding to an age of 1.3 10'° years 
(except in one series where, for comparison purposes, we have 
considered the energy distribution at a burst age of 8 10° years). 
We have explored different values of the ionization parameter 
U, defined as U = yy/cny where yy is the flux of ionizing 
photons impinging a cloud of density ny and c is the speed of 
light. The computations were carried out with the multipurpose 
photoionization-shock code MAPPINGS (Binette et al. 1993: 
BWVMs3; and references therein) which include the effects of 
internal dust not only on the ionization structure but also on the 
emerging line spectrum. We have considered radiation bounded 
models and simply added together the spectra of the two per- 
spectives front and back (‘f’ + ‘b’; cf. BWVMs3) of the ionized 
slab. We did not explore the case of substantial dust opacity 
in the neutral zone of the slab since this would appear an ar- 
bitrary parameter at this stage. We have explored three sets of 


! Demoulin-Ulrich et al. (1984) claimed that the stellar populations 


in elliptical galaxies were unable to produce Ha equivalent widths 
larger than 0.1 A. The reason for such a small value is that they assume 
that the ionizing radiation was produced by horizontal branch stars at 
30000 K (supposed to account for most of the observed radiation at 
1000 A). 


abundances: Z = | (i.e. “solar”), Z = 3 and Z = 1/3, where 
Z represents the total metallicity (dust+gas phase) of the ISM 
with respect with the solar abundance set of Anders & Grevesse 
(1989). The parameter Z is simply a scaling factor of the rel- 
ative abundances of any element (except He) with respect to 
hydrogen. Following the procedure outlined by BWMV3, to 
obtain the gas phase abundance set (Z,,;) actually used in the 
calculations, the total (dust+gas) abundance set was depleted in 
proportion to the dust content of the ionized phase. In all cal- 
culations, we have arbitrarily set the dust content at p = Z/3 
where ju is a parameter describing the dust-to-gas ratio in units 
of the value in the solar neighbourhood. We find that the loss in 
efficiency (relative to recombination Case B) in generating Ha 
due to the internal dust in these particular models at U = 10~‘ is 
only 1.6% and 5% for Z = 1, and 3, respectively. The calculated 
emerging Balmer decrement is 2.79 and 2.83, respectively, and 
is, therefore, not significantly reddened by internal dust. In all 
calculations, we have assumed constant pressure within the slab 
starting with a density of ny = 10° cm~? at the outer face since 
such a density is consistent with the mean value determined 
from the red [S 11] doublet (PJDSB). 


The different panels of Fig.3 display, respectively, the 
computed intensities of [N11] 46583 A, [O1] \6300A, [S11] 
6731 A, [O11] 43727 A, [S m1] A9529 A, [O mt] 45007 A, He tl 
4686 A and [Nett] \3869 A, as a function of U. All the in- 
tensities are given relatively to the Ha intensity. All the dashed 
curves correspond to an age of 1.3 10!° years but for different 
values of Z while the dotted curve corresponds to an age of 
8 10° years (and Z = 1). 


Because photoionization models are better constrained 
when the number of observed lines is larger, we have chosen 
to compare our models with the sample of normal elliptical 
galaxies studied by Demoulin-Ulrich et al. (1984). Indeed, the 
spectrophotometric measurements made by these authors cover 
a wide wavelength range (A\3600-7000 A). We have indicated 
by thick solid lines in Fig. 3 the mean values of the line ratios 
derived from their Table 4 (excluding the radio-galaxy M87). 
No foreground reddening corrections where applied to these ra- 
tios owing to the lack of information on H@. This may well 
imply a substantially higher intrinsic [O 11]/Ha \3727/A6563 A 
ratio. Allowing for this, comparison between the observed and 
computed line ratios consistently indicates a typical value for 
U of 10. 


We now check whether the mean distance to the ionizing 
sources implied by such a value of U is compatible with the 
dimensions of the emitting zones. Assuming that a typical post- 
AGB ionizing star produces 10*’ ionizing photons per second, 
we find that, with a density of 10° cm~?, an ionized cloud must 
lie at a typical distance of 1.6 10'® cm from the star. An inte- 
grated Qy of 10°? photons s~! corresponds to 10° stars. A total 
volume of 10°4/37(1.6 10!8)? cm~3 (~ 6 1075(10? /ny)* kpc?) 
can easily be accommodated within a region of diameter of 
~ Ikpc, even if much flattened or of lower gas density by a 
factor 10. Thus, we consider our model consistent with the ob- 
served excitation level as represented by the line ratios. 
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Fig. 3. Line ratios respective to Ha as a function of the ionization parameter U for models corresponding to photoionization by old post-AGB 
stars. Different curves correspond to different gas abundances (solar, one third solar and 3 times solar), and to different burst ages (8 10” and 
1.3 10'° years). The thick line represents the mean line ratios observed by Demoulin-Ulrich et al. (1984) for the normal ellipticals NGC 2911, 


NGC 2974, NGC 3962, NGC 5077 and NGC 5077A 


The question of the nitrogen overabundance in ellipticals 
can be addressed by comparing the observed [N 11]/Ha ratios 
with the ones computed in our models on the assumption of a 
solar N/O abundance ratio. The study of PJDSB provides a set 
of measured [N !1]/Ha in 95 E and SO galaxies which, except for 
one object, all have [N t!]//Ha < 2.5. The average line ratio is 1.4 
for objects not classified as H i! region-like. Considering that, 
at U = 10-4, our Z = | and Z = 3 models predict an [N 11]/Ha 
ratio of about 1.2, and even the Z = 1/3 model predicts a ratio 
of 0.7, we find little evidence that the N/O ratio in most objects 
should exceed a value of twice solar, at least at this stage of 
our analysis. This seems at variance with the predictions of the 
chemical evolution models of Matteucci & Padovani (1993) for 
elliptical galaxies. However, it would be worthwhile studying 
for normal ellipticals the slit-size effect investigated by Storchi- 
Bergman (1991) for Seyfert 2 galaxies and LINERs. 


5. Conclusion 


We have shown that post-AGB stars in an ageing starburst pro- 
vide about the adequate number of Lyman continuum photons to 
account for the observed emission from ionized gas in elliptical 
galaxies. Assuming a covering factor close to unity, our mod- 
els accounts for both the luminosity and the equivalent width 
of Ha. The observed excitation level of the ionized gas (indi- 
cated by the emission line ratios) can easily be reproduced by 
photoionization models. 


From this study, a general picture emerges, in which ioniza- 
tion of the gas observed in elliptical galaxies would simply be 
provided by their old stellar populations. Our model, of course, 
is prone to revisions, mainly because the evolutionary status of 
low and intermediate mass stars is still unclear. Especially, it 
strongly depends on the stellar mass-loss rates during the var- 
ious evolutionary phases. These mass-loss rates are not well 
known, but certainly, metallicity plays an important role (see 
e.g. a discussion in Brocato et al. 1990). Di Serego Alighieri et 
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al. (1990) have argued that post-AGB evolutionary tracks with- 
out thermal pulse would lead to an even higher production rate of 
ionizing photons than the tracks of Schénberner (1983), which 
were used in the present study. Clearly, much work remains to 
be done before the UV and Lyman continuum output from old 
stellar populations is understood. 

Within the framework of our model, further studies are also 
needed, and in particular, detailed observations and modelling 
of individual objects in order to reproduce simultaneously the 
observed continuum — including measurements shortward of 
2000 A — and the emission lines. It is also important to check 
whether additional sources of ionization are required in specific 
cases, or, on the contrary, whether part of this ionizing radiation 
may escape into the intergalactic medium. 

We find no evidence for an N/O abundance ratio exceeding 
a factor two that of solar. Although the observational constraints 
are still scarce, it seems that the chemical composition of the 
ionized gas is close to solar (within a factor 3). 

Our model also partly applies to LINERs, a phenomenon 
which is stronger and better studied in later type galaxies (Keel 
1983a,b). We are not proposing that all the emission in LIN- 
ERs can be explained in this way but simply that there exists a 
minimum level of ionizing photons produced in every evolved 
stellar system, which has to be taken into account. An active 
nucleus in some ellipticals is not ruled out, however, especially 
in order to power the radio emission. 
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Abstract. Bidimensional spectroscopy of the distant 
(z = 0.471) radio-galaxy 3C 435A has been performed with 
the integral field spectrograph TIGER at the CFHT. At least 
five components are identified around the central one. The 
[OI] 3727A and [O III] 4959,5007A maps are remarkably su- 
perimposed with the 1.4 GHz radio lobe isophotes, showing the 
narrow relationship between the radio plasma and the ionised 
component, while the 4.8 GHz radio emission traces a shock 
region. The [OII]/[O III] ratio map follows the external part of 
the radio lobe rather than the central jet, and is strongly distorted 
by multiple components. This high number density is confirmed 
by the velocity field. The maximum ratio value (~ 8) indicates 
a low level of ionisation. For the first time, stellar and nebu- 
lar components are separated at a large redshift, indicating that 
the alignment effect is mainly due to nebular [O II] and [O IIT] 
emission. With the help of our spectrophotometric model, the 
stellar continuum of the central component A is well fitted with 
an old (more than 15 Gyrs at z = 0) elliptical galaxy model, 
while nebular lines witness a faint star formation process. Two 
companions B and D, located near the radioaxis, fit evolved 
(= 10 Gyrs) models, ruling out a current star formation activ- 
ity induced by the radio-jet. Another companion, C, located in 
a direction roughly perpendicular to the radioaxis, has a blue 
color of irregular type, a few Gyr age, detectable emission lines 
and high velocity. Component E, roughly aligned with A and C, 
is essentially emitting in the [OII] line. Masses and magnitudes 
of some components are given and, among various radiogalaxy 
models, a scenario of ionised gas cocoon due to overpressure, 
or to interactions of the relativistic electrons with the intergalac- 
tic medium, with a possible axis-perpendicular gaseous disk, is 
favored. 


Key words: galaxies: 3C 435A; evolution; structure; jets — radio 
continuum: galaxies 
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* based on observations collected at Canada-France-Hawaii Tele- 
scope, which is operated by CNRS of France, NRC of Canada, and the 
University of Hawaii 


1. Introduction 


Radiogalaxies are the most distant objects in which stellar popu- 
lation signatures can be detected. So they are preferential targets 
for understanding primeval star formation in galaxies, providing 
constraints on theoretical galaxy formation models. 


The complexity of radiogalaxy structures makes difficult 
the detection of star formation signatures, as well as the sep- 
aration of nebular and stellar components at low spatial reso- 
lution. The radiogalaxy B2 0902+34 at z = 3.395 (Lilly 1988) 
is a significant example; the important gap from the far-UV 
to visible (rest frame) was thought to originate from stellar 
emission, leading to an age greater than | Gyr. In fact, huge 
{OIII] and HG emission lines, discovered with the CGS4 in- 
strument at the UKIRT telescope (Eales et al. 1993), are bright 
enough to explain this gap by a nebular rather than stellar ori- 
gin. If the emitting gas is photoionised, this may be a sign of 
intense star formation, and could indicate a very young galaxy 
(age lower than 0.5 Gyr), while stellar continua are tracers 
of past star formation (Rocca-Volmerange 1993). In fact, the 
origin of emission lines and their evolution with redshift are 
still debated. Many scenarii to power emission lines are pro- 
posed, involving photoionisation, shocks or scattering. More- 
over, it is not clear how emission lines or continuum could 
explain the optical and radio axis alignement, discovered in 
a significant sample of radio-galaxies (McCarthy et al. 1987; 
Chambers et al. 1987). Models successively attributed such 
an alignment effect to young stars, interactions between 
relativistic electrons and the ambiant plasma, scattering 
by dust, hot or cold electrons, or inverse Compton scat- 
tering (Rees 1989; de Young 1989; Begelman & Cioffi 1989; 
Daly 1992a; Daly 1992b; Bithell & Rees 1990). Bidimensional 
spectroscopy appears to be the best way for identifying the basic 
physical processes driving star formation in distant radiogalax- 
ies. Another constraint is given by the low level of polarisation 
measured by di Serego Alighieri et al. 1993. 


The distant (z = 0.471) radiogalaxy 3C 435A has been ob- 
served with the integral field spectrograph TIGER at the 3.60 m 
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Fig. 1. Direct R image of 3C 435A, with isophotes 68% to 100%, step 


2%. Component A is at (0,0); locations of the six other components (B 
to G) are labeled. Coordinates are arcseconds 





Canada France Hawaii Telescope. Observations and data pro- 
cessing are presented in Sect. 2, while results and discussion 
are given in Sect. 3 and 4. 


2. Observations and data reduction 


The radiosource 3C 435A, identified with an R~19 galaxy at 
z=0.471 (Smith & Spinrad 1980), belongs to a projection pair 
with 3C 435B at z=0.865. It was selected for its rich clus- 
ter environment (Yates et al. 1989) and its redshift allowing to 
observe both [O II] 3727A, HG and [OIII] 4959,5007A emis- 
sion lines. The observations were carried out in July 1992 with 
the TIGER integral field spectrograph (Courtés G. et al. 1987; 
Pécontal 1991), installed at the 3.60m CFHT, with a mean 
0.7” seeing. The spatial sampling is performed by a micro-lens 
array. Each spatial element then gives a spectrum on the CCD, 
through a classical spectrograph. Spectra are extracted from the 
CCD image using a dedicated software, developed at Obser- 
vatoire de Lyon, and integrated into the ESO MIDAS package 
(Rousset 1992). 

In direct imaging mode, the spatial sampling was 0.16”, 
and the resolution seeing-limited over the 38” x36” field. Af- 
ter a light median filtering, an effective resolution of 1’ was 
measured. Seven direct imaging frames correspond to a total 
integration time of 1620 s. The astrometric accuracy, from the 
point-like nucleus of 3C 435B (Smith & Spinrad 1980) and an 
”A” star (Wills B. & Wills D. private communication) is + 0.4”. 

Each spectrographic frame covers an 11” x11” field, with 
a 0.61” spatial sampling, giving about 350 spectra with a 8A 
spectral sampling. Observations were made within two wave- 
length ranges: 5000A-7000A (hereafter called ”V” configura- 
tion), and 6500A-8500A (hereafter called ”’R” configuration) 
with two exposures in V (2x3600s), four in R (2x 3600 s, 
2463 s and 1808 s). 


Data reduction was performed using the TIGER software 
package. Wavelength calibration was checked with the [OT] 
5577.4A night sky line, yielding a + 1 A measured accuracy. 
Using the same line, we measured an effective final spectral 
resolution of 20.8 A. 

Data have been spatially smoothed (averaged over neigh- 
bouring lenses) to increase the S/N ratio of spectra. Fitting of 
the [OII] and [OTII] emission lines with Gaussian profiles on 
each spectrum gave us wavelengths and continuum-free inten- 
sities; this allowed us to build intensity, line ratio, and velocity 
maps. 


3. Results 
3.1. Morphology 


The large scale structure, observed on a direct short ex- 
posure image (Fig. 1), confirms the existence of the pre- 
viously known 3C435A components (McCarthy et al. 1989; 
di Serego Alighieri et al. 1993), while several additional ones 
are identified hereafter, owing to our 2D spectrographic data. 
At a smaller scale, images are rebuilt by integrating spectrum 
continua over broad-band ”R” and ”V” wavelength intervals 
(which are not the standard Johnson filters) and mapping [O IT] 
and [O III] energies (from fitted gaussian areas). Components A 
and B show smooth continuum isophotes in both V and R, while 
components C and D are too faint to show any structure. New 
nebular components F and possibly E clearly appear through the 
[O Il] and [OTI)/[O III] ratio maps (Fig. 2). The faintest [O IT] 
isophotes follow the crescent structure of the 1.4 GHz radio 
lobe (Fig. 5). The [O III] emission is largely extended towards 
North in the hot spot region and, at a lower level, towards West 
in direction of component C. The [O II]/[O III] ratio map shows 
that the central isophotes are distorted by the D and G compo- 
nents (Fig. 2), confirmed below by the velocity field; the line 
ratio peaks at ~ 8 along a crescent isophote following the radio 
structure. 


3.2. Spectroscopy 


The six components A, B, C, E, F and G show stellar and/or 
nebular lines which all confirm their membership to the radio- 
galaxy. Only two of them, A and B, have a S/N ratio high enough 
for performing a comparison with evolutionary models. Table 1 
gives apparent V and R fluxes and [OII] and [O II] emissivi- 
ties for each component; uncertainties are essentially due to flat 
field correction. 


— Component A (Fig. 4) has a significant continuous stellar 
energy distribution from 5400 A to 8500 A; a fairly large 
(900 A) overlap between V and R gives a check of the data 
processing validity. A precise estimate of the redshift gives 
z = 0.470 + 0.003 from both [O II] and [OTI]. Apparent 
magnitudes are V=21.5 and R=20 + 0.2, without reddening 
correction. 
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Fig. 2. Published radio maps, and maps rebuilt from TIGER spectroscopic data. Coordinates are arcseconds, and components A to G are labeled. 
Reconstructed map isophotes are in 10~'* erg s~' cm~? arcsec? except for the line ratio map (unitless). Radio maps: contours in mJy; see 
text for details. R-band continuum map: continuum flux in a 7500/1900 A band. Plotted isophotes: 5, 10 and 15 to 95 with a step of 10. [OI] 
3727 A map: line flux derived from a gaussian fit. Plotted isophotes: 0.01 (dashed line) and 0.05 to 0.40 with a 0.05 step. [O III] 5007 A map: 
line flux derived from a gaussian fit. Plotted isophotes: 0.01 (dashed line) and 0.10 to 1.70 with a 0.20 step. [O II] 3727 A/{O111]5007 A map: 
Line ratio map. Plotted isophotes: 2.5 to 8.5 with a step of 1 


— Component B (Fig. 4) has magnitudes of V=22.3 and R= 
20.9 + 0.2 and shows an emission-free continuum. It is cen- 
tred on -1.8”, 4.7”, A being at (0,0). 

— Component C shows faint S/N ratio continuum and emission 
lines shifted by 600 kms~' relative to the central compo- 
nent. About magnitudes, a rough value of 23, in both V 


and R, confirms the irregular galaxy nature predicted by the 
noisy spectrum. 

— Component D , too faint in the V band, is marginally de- 
tected in R, with a magnitude of 24.5 within a 1.8” circle 
centred on position (-2.3”,1.6”). 

— Component E is very close to A; if the [O IT] emission tightly 
follows the radio lobe it could be just a projection effect 
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Table 1. V (6000/2000 A) and R (7500/1900 A) continua fluxes, and integrated [OII] 3727 A and [OIII] 5007 A line fluxes. Units are 
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Components V flux R flux [OIT] flux 


Accuracy 


[OTT] flux 


Other lines 





A+E 1.8 
0.9 


0.5 


4.0 
1.4 
0.2 


~ 


20% 
30% 
50% 


0.21 0.60 





Call(H+K), Mgb,MgH,HB 
NeV, Nelll, HB, Hy 





due to the curvature of the lobe. The [OII] line emission 
maximum is centred at (-1.6”, -0.2’"), while [O III] is not 
detected and [Hel] is possible. Its velocity is estimated to 
280 kms~' relative to A (e.g. z = 0.469). 

— Component F, centred at (3.2” 3.4”), is located in the inter- 
action zone of the radio jet with the intergalactic medium. 
The velocity is estimated to 180 kms~!. The high spectral 
index of the UV continuum traces a non-thermal emission. 

— Possible component G merely reaches the detection thresh- 
old, but will be confirmed below by the velocity field. 


3.3. The velocity field 


Strong negative velocities, relative to the A central component 
(Fig. 3) are witnesses of gas efficiently blown out by the radio 
jet, and suddenly stopped in the hot spot zone F by the intergalac- 
tic medium. As shown below, [O II] isophotes follow external ra- 
dio isophotes; there is no intensity enhancement along the radio 
jet, but only a patchy structure due to the presence of components 
D and G. These two negative velocity components may witness 
the presence of clouds or inhomogeneities located between the 
radio plasma and the ambiant gas, which do not affect the radio 
emission but compress the ionised gas. This is similarly found in 
3C326.1, another distant radiogalaxy (McCarthy et al. 1987b). 
An interesting point is the quite symmetrical value of pos- 
itive and negative velocities (~ +400 kms~ !) which corre- 
sponds to a total Av~800km.s~! for a global size of 100kpc 
(Ho=50km.s~'.Mpe~!, qo=0.5, Ao=0), assuming a symmetri- 
cal southern lobe. 


3.4. Comparison with radio and polarimetric data 


A comparison of the [OII] and [OTII] distributions with 
the radio maps is shown in Fig. 5, respectively at 1.4 and 
4.8 GHz. Nebular isophotes follow the radio plasma along 
the radio lobe, giving interesting constraints on the physical 
origin of the light emission. The comparison of the [OTII] 
map with the 4.8 GHz map (J. Hutchings, private communi- 
cation) shows that the NW emission of the [OIII] map is at 
the edge of a compact radio component. The A-F distance 
gives a size of 30kpe (Hp=50km.s~!.Mpc™', qo=0.5, Ao=0). 
The southern part of the radiogalaxy has not yet been ob- 
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Fig. 3. [OII] velocity field v({O II]) - vo (where vo is the velocity of 


component A) The total range is [-380,410], with a 80 km.s~! step. 
Dashed isovelocities are for negative values 





served with TIGER. The polarisation measurements published 
for 3C 435A (di Serego Alighieri et al. 1993) show a polarisa- 
tion factor lower than 5.7% for B and lower than 3.6% for all the 
components together. The electric vector, with a position angle 
of 138+15° for the object as a whole, is roughly perpendicu- 
lar to the radio axis (Cimatti et al. 1994). This low polarisation 
could be related to a low [OII] emission along the radio-jet. 


4. Discussion 


Among the most pressing issues concerning distant radiogalax- 
ies, the two questions of confinement and age will be decisive 
for understanding their formation and evolution. From our anal- 
ysis, nebular and stellar emissions as well as velocities of the 
ionised [OI] gas are identified. Moreover the [O II]/[O III] ra- 
tio gives clues on the nature of nebular emission while galactic 
evolution models determine ages of stellar populations. 


4.1. The alignement effect: a nebular origin 


The alignment effect was discovered, on a large sample of dis- 
tant galaxies, to increase with redshift (McCarthy et al. 1987; 
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Fig. 4. Fits of components A and B spectra (shifted to z = 0) with the Synthetic Atlas of Galaxies. Abscissa: A . Ordinates: arbitrary scales in 
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Chambers et al. 1987). For this reason, it was claimed to origi- 
nate from physical processes distinct from those involved at low 
redshift (van Breugel et al. 1985; Baum & Heckman 1989). 
Some models interpret this alignment by a star formation 
activity (Rees 1989; de Young 1989; Begelman & Cioffi 1989; 
Daly 1990; Daly 1992a; Bithell & Rees 1990). Star formation 
was proposed to occur in cold clouds inside radio lobes, where 
compression induces gravitational instability, triggering star for- 
mation (Rees 1989; Bithell & Rees 1990). From our spectral 
data, this scenario could only be possible into D and G com- 
ponents, but the faint spectrum of D is red and corresponds to 
an evolved population, while in G no continuum reaches the 
detection threshold. Another possibility is the scattering of the 
anisotropic radiation emerging from the galaxy nucleus, induced 
by dust, hot or cold electrons (see Daly 1992a; R&ttgering 1992 
for a review). More recently (Daly 1992a; Daly 1992b), an in- 
verse Compton scattering of microwave background photons 
by relatively cool relativistic electrons has been proposed. In 


these two last cases, the upscattered radiation would pref- 
erentially be aligned with the radio axis. For 3C 435A, the 
optical axis was previously defined by components A and 
B (di Serego Alighieri et al. 1993) and was therefore poorly 
aligned (within 25°) with the radio axis. In fact, a significant 
alignment is visible in [O II] between A and F and corresponds, 
within 5°, to the 1.4GHz radio axis. Analysing the [OII] and 
[O III] maps, one sees no evidence of scattering along the ra- 
dio axis. [OIII] is maximum in the central component A and 
in F, tracing a zone of strong interaction with the intergalactic 
medium, as shown by the 4.8 GHz radio map. The [O IT]/[O I] 
ratio is more significant because its maximum value (~8) fol- 
lows a crescent distribution along the radio isophotes, confirm- 
ing the narrow relationship between radio and ionised gas. This 
ratio indicates a gas density of 10° to 10* cm~3, with a tem- 
perature of 4000 K and a low ionisation parameter of 10~* 
(Osterbrock 1974; Stazifiska 1984). This distribution, projected 
on the line of sight, would agree with the overpressured cocoon 
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Fig. 5. Left: gray-scale [O II] 3727 A reconstructed map with the isophotes of the 1.4 GHz VLA map of McCarthy et al. 1989 overimposed. Right: 
gray-scale [O III] 5007 A reconstructed map with the isophotes of the 4.8 GHz VLA map (Hutchings J., private communication) overimposed. 
Strange features in the radio isophotes are only artefacts of the scanning process 


model (Begelman & Cioffi 1989), or with interaction between 
relativistic electrons and ambiant gas. 


4.2. Overpressured cocoon and ionised disk 


This cocoon model interprets the [O IT] and [O III] emitting ma- 
terial as a high pressure zone of shocked gas produced by the in- 
teraction of the radio plasma with the intergalactic medium. The 
pressure plays a fundamental role, confining the radio plasma 
without the need of any magnetic field. The observed structure 
(Fig. 5), with an [OIII] emission zone in front of the 4.8 GHz 
hot spot at the end of the radio jet, and peripheral [O II] emis- 
sion closely following the 1.4 GHz radio isophotes, agrees with 
a bow shock hypothesis. However such a model does not fit 
the radial expansion near the radiogalaxy itself. Observations 
of component A surrounded by high velocity extensions would 
be in better agreement with a disk of [O II] emitting ionised gas, 
roughly perpendicular to the radiojet. Such an image is remi- 
niscent of the large scale starburst structure observed in M82 
(Tully, private communication). 

The [OII] ionised cocoon interpretation is also in agreement 
with the unified model of powerful radiogalaxies and radio-loud 
quasars (Barthel 1989). The [OII] external distribution would 
explain by projection the apparent isotropy for the two classes 
of objects (Hes et al. 1993). A more precise modelling is in 
progress. 


4.3. Ages of the various components 


According to the complex structure of a radiogalaxy, we 
have to define the age of various stellar components, the 
non-thermal process being much more difficult to date, 
even if its lifetime is shorter. Fits of components A and 
B (Fig. 4) were computed using our Atlas of Synthetic 
Galaxies (Rocca-Volmerange et al. 1988) and a x? minimum 
test. Components A and B have respective ages of 11 and 


8.5 Gyrs at z=0.47; this would mean respective ages of 17 
and 14.5 Gyrs at z=0, for extreme cosmological parameters 
(Ho=50km.s~!.Mpc™! and qo=0). That corresponds to the stel- 
lar population of globular clusters. Clearly component A is the 
oldest. A flat Universe is not acceptable. 

If we adopt a MIL ratio of 10 at 10 Gyrs for our ellipti- 
cal model, this implies M4=4.6 10'° Mo. For the B compo- 
nent, we have M/L=8 at 8.5 Gyrs with Mp=1.1 10'° Mo. At 
the opposite, component C, which is barely detected, is much 
younger, no more than 3Gyrs at z=0.47. A MIL ratio of 5 gives 
Mc=0.85 10'° Mo. Large ages for A, B and D components lo- 
cated inside the radio axis cone, and small age for component 
C, suggest constraints for galaxy formation models; star forma- 
tion was active along the radio jet at the earliest epoch while the 
nucleus activity is much more recent, and is not triggering a sig- 
nificative star formation. A possible consequence is a recurrent 
process for active nuclei (Haenhelt & Rees 1993). 


4.4. The high number density of neighbours 


A high number density of components is observed in the nearby 
environment. The central component A, companions B and 
(likely) D are old galaxies, while C has a young population. 
G, detected in [OII] nebular emission, is below the detection 
threshold in stellar light. Component E, identified by a strong 
[O II] emission, has a high velocity relative to A, and is located 
less than 1” from it. This could be typical of an interaction pro- 
cess, its gas fueling the active nucleus A. We also consider the 
possibility to observe the projection effect of a nebular cocoon, 
or disk (see Sect. 4.2) Component F is essentially a nebular one, 
where strong [OIII] emission could be triggered by shocks. 


5. Conclusion 


Such two-dimension observations show that the alignment 
effect is essentially due to the nebular emission (as in 





26 B. Rocca-Volmerange et al.: Structure and age of 3C 435A 


van Breugel et al. 1985). [OII] isophotes, following external 
part of the radiolobe, favor the scenario of gas being overpres- 
sured by radio plasma expansion, or the interaction of relativis- 
tic electrons with the ambiant gas. This could partly explain the 
confinment of the radio plasma. 

There is no evidence of current star formation from emission 
lines except in component A. Datation of spectra gives large 
stellar population ages in A and B, similar to globular cluster 
ones. That rules out models of star formation triggered by the 
radiojet. A high object number density is observed in the nearby 
environment, implying distorted isophotes for the radiogalaxy 
3C 435A. A larger sample is needed for improving constraints 
on galaxy formation models. 
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Abstract. The intracluster gas of clusters of galaxies shows a 
metal enrichment with nearly solar abundances of metals. Two 
possible mechanisms are indicated for the injection of material 
by member galaxies: ram pressure stripping and protogalac- 
tic winds. In this article it is proposed a new static model that 
with respect to the non-polytropic model formulated by Cav- 
aliere & Fusco-Femiano (1981), utilized to jointly fit all the 
available X-ray and optical imaging and spectroscopic data of 
the Coma cluster (Fusco-Femiano & Hughes 1993), includes 
also the contribution of supernovae-driven protogalactic winds 
to the heating and mass of the intracluster medium. The use of 
this model for the analysis of spatially resolved X-ray spectra 
should allow to discriminate between the wind and ram-pressure 
stripping scenario, regarding the injected energy and the metal 
enrichment of the intracluster gas. 


Key words: galaxies: clusters of — intergalactic medium — galax- 
ies: evolution 





1. Introduction 


Fusco-Femiano & Hughes (1993; hereafter FFH) have shown 
that the non-polytropic model (Cavaliere & Fusco-Femiano 
1981; hereafter CFF), based on the assumption that the intr- 
acluster gas temperature is proportional to the square of the 
galaxy velocity dispersion everywhere throughout the cluster, 
is fully consistent with the density and temperature distribu- 
tions of the gas and galaxies in the Coma cluster. From the joint 
analysis of the X-ray and optical data, FFH were able to de- 
rive all the parameters that characterize the model. This model 
has been taken in consideration after that Hughes et al. (1988) 
pointed out that self-consistent polytropic models were unable 
to provide acceptable fits to the X-ray imaging and spectral data 
on Coma cluster. The validity of the CFF non-polytropic model 
based on the assumption that kT(R) ~ o*(R) suggests that a 
significant fraction of the thermal energy contained in the intra- 
cluster medium (ICM) of the Coma cluster (as well as other rich 
galaxy clusters) may have come from the interaction between 
the galaxies and the ambient cluster medium. 
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White (1991) has studied the energetics associated with 
metal enrichment of the intracluster gas. Galaxies may inject 
metals into ICM via winds or by the stripping process. He shows 
that while the energy input from gas stripping process enters as 
o”, energy input from galactic winds depends on the depth of the 
protogalactic well of the individual galaxy as well as the wind 
speed, and not only on o”. If the injection is due to supernovae- 
driven protogalactic winds, the observed abundances reflect the 
degree of stellar processing because the metals are distributed 
fairly homogeneously. Otherwise, in the case of centrally con- 
centrated stripping, the overall abundances could be much less 
than those observed in cluster cores. 

A related problem regards the ”G-discrepancy” of X-ray 
clusters where 3 = umyo?/kT (u ~ 0.6 is the mean molecu- 
lar weight and my is the proton mass) is a measure of the ratio 
of energy per unit mass in galaxies to that in gas. The hydro- 
static isothermal G-model (Cavaliere & Fusco-Femiano 1976) 
has been used successfully to fit the X-ray surface brightness 
emission from numerous clusters of galaxies; a value (3; ;z is de- 
rived (Gorenstein et al. 1978; Branduardi-Raymont et al. 1981; 
Abramopoulos & Ku 1983; Jones & Forman 1984). Jones & 
Forman (1984) find from Einstein X-ray surface brightness data 
for clusters an average value of (Git) ~ 0.64 indicating that 
the ICM is hotter and more extended than galaxies. Instead, di- 
rect spectroscopic measurements of the gas temperatures and 
optically determined galaxy velocity dispersion produce val- 
ues Bsnec ~ 1.2 (Mushotzky 1984, 1988). The interpretation of 
this value is the opposite of the previous one (';;z. Evrard (1990) 
finds the same discrepancy performing a hydrodynamic simula- 
tion of the ICM. This is attributed to incomplete thermalization 
of the gas in the ICM. From the analyses of the EXOSAT cluster 
observations, Edge & Stewart (1991) report that with the im- 
provement in the optical data, ¢;¢ and Gspe- are both less than 
unity resolving in this way the ”G — problem”. The better qual- 
ity of the optical data has reduced the high velocity dispersions 
in several clusters due to the presence of significant velocity 
substructures (Fitchett & Smail 1991) and {s,e-. > 1 may be 
an indicator of the presence of such substructures in clusters of 
galaxies. 

Open problems remain the value of 6 < 1, that involves 
an additional energy source for the gas respect to the galaxies, 
and an increasing value of @ with cluster mass (White 1991; 
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David et al. 1991); the best fit values of G range from ~ 0.41 
for groups (Kriss et al. 1983) to G ~ 0.64 for rich clusters (Jones 
& Forman 1984). 

In this paper we introduce in the CFF non-polytropic model 
the possibility that protogalactic winds may contribution to the 
heating of the ICM along with the injection of metals. In this 
model we consider that the major sources of intracluster gas are 
the primordial gas, the gas stripped from the galaxies and the 
gas injected by galactic winds. 

The outline of the paper is as follows. In Sect. 2 we discuss 
the features of the non-polytropic model with galactic winds. 
Sect. 3 reports the implications that this model has on the fits to 
the brightness distributions of poor and rich clusters with respect 
to those obtained using the isothermal 3-model or the CFF non- 
polytropic model. Discussion on the meaning of 6 < | and of 
the increase of 3+;, from poor to rich clusters appears in Sect. 
4 and conclusions are in Sect. 5. 


2. The model 


The CFF non-polytropic model assumes that the intracluster 
medium is composite of primordial gas that is compressed and 
heated by gravity during the cluster collapse and by gas stripped 
from galaxies by the ambient medium. The kinetic energy of the 
ejected gas would be degraded into internal energy by collisions 
with the ambient intracluster gas. Both these two components 
of the ICM have a specific energy roughly equal to 3/207. Most 
of the ICM must be primordial since the gas mass is much 
greater than the observed stellar mass of the cluster galaxies 
(Blumenthal et al. 1984; David et al. 1990). In this model the 
relation 


kT(R) = pmyo?(R)/B (1) 


holds everywhere throughout the cluster. Two objections at this 
relation are that the heating of the primordial gas is a very com- 
plex issue (the processes that heat this gas would be probably 
different from those which would have heated the galaxies) and 
that simulations of cluster formation indicate that the galaxy ve- 
locity anisotropy should vary with radius. However, FFH have 
demonstrated, for the first time, how a physically- motivated 
static model in which the heating is controlled by local condi- 
tions, can jointly fit all the actual X-ray and optical data on the 
Coma cluster. 

Inserting the relation (1) into the equation that governs the 
equilibrium condition for the gas and the galaxies in the same 
potential well 


umyn/ pg = dp/dpg (2) 


it is possible to obtain the distributions for the gas density and 
temperature taking into account that p = kT'n and pg = pgo? 
(Cavaliere & Fusco-Femiano 1978; CFF): 
n/N = (pa/Pao)*(o?/02)>—| (3a) 


T/T. = 07/02 = Is/2(W)/Is2(Wo) * I3/2(Wo)/Is/2(W) (36) 


where 
w 
I,(W) = [ e "n%dn 
0 


is the incomplete gamma function +(q + 1, W). 
The galaxy density is given by 


pa/ Pao = eW—W* Ix(W)/I5/2(Wo) 


In turn the normalized gravitational potential 


Ri 
win=G | M(r)r~?0~?dr 
R 


is derived from the full Poisson equation that includes the total 
density p(W) = n(W) + pg(W) (any "virial mass” distributed 
like the galaxies is included in pg (King 1966; Cavaliere & 
Fusco-Femiano 1978)). The cluster boundary R; is defined to 
be the radius where pg = o = 0. 

The enrichment with heavy metals of the intracluster gas 
may be obtained by stripping or by galactic winds. So, a mass 
fraction of the ICM and an higher gas temperature can derive 
from ejection by galactic winds. White (1991) has shown that 
the energy input from galactic winds depends on the depth of 
the potential well of the individual galaxy as well as the wind 
speed, and not only on o°. 

With respect to the CFF non-polytropic model, we introduce 
in this model the possibility that a contribution to the heating of 
the ICM may derive from the ejection of material via galactic 
winds. The relation (1) is modified including this effect 


kT(R) = pmy(o?(R) + Ewind)/B 
where €,,inq takes into account the heating by galactic winds. 


Inserting this new relation in the equation (2), the gas density 
and temperature distributions are now the following: 


n/No = (pa/ Pao) [(o?/o2 + q)/(1 + qo 'e~ fun (4a) 


T/T = (07/02 + q)/(1+q) 


where q = €wina/o? and 


(46) 


W (1 +e-*29/?/Ty9(x)) dx 


W. (0? /02)+4q 


fwina = 98 





These gas density and temperature distributions derive from the 
contributions to the ICM of the primordial gas, of the stripped 
gas from the galaxies and of the material ejected via protogalac- 
tic winds. In this model the relevant parameters, which are to be 
determined by fits to the X-ray and optical data, are the follow- 
ing: No, the central gas density; T,, the central gas temperature; 
R,, the core radius of the galaxy space distribution; W,, the 
central value of the normalized gravitational potential; G, the 
ratio of the specific energy between the galaxies and gas; and 
€wind the energy deposited by the galactic winds in the ICM. 
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The expression that allows computation of the cluster bind- 
ing mass is derived by combining the spherically symmetric 
condition of hydrostatic equilibrium with the ideal gas law: 


kT (2 ‘ =) __GM(R) 


umy \n TT) R2 -" 


which may be written as 


dn | dT \ dW 
ndw Taw) © “dy 





M(y) = — 
(y) mal 


kT ( 
where y = R/R,. For this model the derivatives assume the 
following forms 


dn eWwis2 e—Wwy3/2 
=e, = 6 +(6 ~ I) + - 
ndWw Ts/2(W) T3/2(W) 
dT —_whp(W)W%!? — Is(W)W?/? 


TdW Is(W)Is/2(W) 


dfwind 
dw 








where dfwing/dW = qBF(W) with 


1+e-WW?/?/T,.(W) 
o*/o3+4q 





F(W)= 


The gravitational gradient dW//dy is given by the solution of 
the full Poisson equation (Cavaliere & Fusco-Femiano 1978). 

In the next section we discuss the behaviour of this model 
regarding the fits to the brightness distributions of poor and rich 
clusters of galaxies. 


3. Characteristics of the model 


The analytical expression given by 


S(b) = S(O) + (6/ReYP PO? (5) 
derived from the hydrostatic isothermal G-model (Cavaliere & 
Fusco-Femiano 1976), has been used succesfully to fit the X-ray 
surface brightness emission from numerous clusters of galax- 
ies (Gorenstein et al. 1978; Branduardi-Raymont et al. 1981; 
Abramopoulos & Ku 1983; Jones & Forman 1984). In this re- 
lation b is the projected radius and R, is the core radius. This 
model represents the limit, y = 1, of the polytropic models. 
The fits to the X-ray surface brightness distributions of clus- 
ters allow the determination of the gas density, while informa- 
tions on the temperature distributions has been more difficult to 
obtain due to technical limitations of instrumentation flown on 
past X-ray Astronomy satellites. The Coma cluster is one of the 
few clusters for which there is some knowledge of the temper- 
ature distribution. Watt et al. (1992) has shown that the Coma 
cluster is nearly isothermal out to roughly 30’, while X-ray spec- 
tral data obtained from collimated instruments with field of view 
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Fig. la. Temperature distribution derived from the fit to all the avail- 
able X-ray and optical data on the Coma cluster using the CFF non- 
polytropic model (FFH) 
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Fig. 1b. Density distribution derived from the fit to all the available X- 
ray and optical data on the Coma cluster using the CFF non-polytropic 
model 


from 45’ to 3°, strongly require a decrease of the temperature 
in the cluster atmosphere. 

The CFF non-polytropic model was used to fit all the avail- 
able X-ray data on Coma cluster (FFH), after that Hughes 
et al. (1988) pointed out that self- consistent polytropic mod- 
els were unable to provide acceptable fits to the X-ray imag- 
ing and spectral data on Coma. From the joint analysis of the 
X-ray data (namely, X-ray surface brightness from the IPC of 
Einstein , spectral data obtained by the X-ray satellites Tenma, 
EXOSAT and Ginga ) with the optical data (namely, surface 
galaxy density and profile of the galaxy velocity dispersion) 
FFH derived the distributions of the gas density and tempera- 
ture with the following parameter values : 6 = (0.59 — 0.68); 
R, = (7'.4 — 9'.2); W, = (8.2 — 9.3); T, = (8.8 — 9.4)keV 
and n, = (2.57 — 3.03) x 10~3cm~3. The central velocity 
dispersion, o,, is a derived parameter based on the definition 
of (; it lies in the interval (900-1000)km/s. Figure 1a,b shows, 
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Fig. 2. Relations between the parameter q = €wina/o2 and G for differ- 
ent values of @+i¢ derived using models with q = 0. 3¢;+=0.41,0.64,0.9 
for curves a,b,c, respectively 
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Fig. 3. This relation shows how to find the angular coefficients of the 
straight lines of Fig. 2 for different values of 3+;, derived using models 
with q =0 


respectively, the temperature and density distributions of the in- 
tracluster gas for parameter values taken at the midpoint of the 
above intervals. 

FFH demonstrated that the density profile of the CFF non- 
polytropic model that gives an excellent fit to the surface bright- 
ness distribution of Coma is virtually indistinguishable from that 
of the isothermal G-model n/n, = (pq/pao)? due to the slow 
radial decrease of the galaxy velocity dispersion in equation 
(3a). The excellent quality of the fit obtained by FFH with the 
CFF non-polytropic model shows clearly that the spatial data 
can be fitted satisfactorily without requiring that the gas in the 
central region of the Coma cluster be in a fully isothermal state, 
as suggested by previous analyses with polytropic models (in 
particular, Jones & Forman 1984 and Hughes et al. 1988 used 
the isothermal (-model, see Eq. (5)). It is relevant to remember 
that because of the limited energy band width of the IPC, the 
surface brightness data are quite insensitive to the value of the 
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Fig. 4. Temperature,(1), and density,(2), distributions normalized to 
their central values for different models: a) isothermal G-model (G = 
0.64); b) CFF non-polytropic model (G = 0.64;W, = 8.75; R. = 
8'.3;T, = 9keV;n. = 2.77 x 10~%cm™*); c) This model (8 = 
0.83;q = 0.3; the values of the remaining parameters are the same 
of previous model b)) 


central temperature J, that, in the spatial analysis of FFH, was 
set to a value of 9 KeV and kept it fixed. 

The introduction of the contribution of galactic winds to the 
ICM in terms of injection of material and heating, involves an 
increase of the gas density and temperature. In order to have 
the same distribution of the gas density, derived from the CFF 
non-polytropic model by FFH (Fig. 1b), able to reproduce the 
surface brightness data of the Coma cluster, is necessary to 
increase the value of @. Figure 2 shows the relation between 
g= €wina/o2 and ( that must be interpreted in the following 
way: the contribution of the protogalactic winds to the energetic 
and mass fraction of the ICM involves an increase of (jz in 
order to have the same density distribution derived from previ- 
ous models (isothermal G-model or CFF non-polytropic model). 
The increase of (+, with respect to the values derived with 
models in which gq = 0 depends linearly from the value of the 
quantity q. Figure 3 reports the relation that allows to find the an- 
gular coefficient of the straight lines of Fig. 2, for a given value 
of 3¢iz, obtained using the isothermal 6-model or the CFF non- 
polytropic model (q = 0). Figure 4 shows that this model with 
protogalactic winds that reproduces the same density distribu- 
tion of the isothermal G-model or of the CFF non-polytropic 
model, able to fit the brightness distribution of the Coma cluster 
(Hughes et al. 1988; FFH), has a different temperature distribu- 
tion with respect to both these models. 

The galactic winds heat the intracluster medium to temper- 
ature greater than that associated with the velocity dispersion of 
the galaxies, namely the temperature produced by gravitational 
heating and ram stripping. This different temperature distribu- 
tion involves a different binding mass distribution, as shown in 
Fig. 5, with respect to the models in which q=0. 

The increase of the value of 6;¢ and of the binding mass 
depends on the energy injected in the intracluster medium by the 
galactic winds. In order to estimate this energy it is necessary to 





R. Fusco-Femiano: Non-Polytropic model with galactic winds for clusters of galaxies 31 





MASS/10*#15 (solor masses) 








0.00 T T T T T T T T 
.00 = .50 





T 

1.00 1.50 2.00 2.50 3.00 3.50 4.00 4.50 5.00 
R (Mpc) 

Fig. 5. Binding mass (1) and gas mass (2). a) CFF non-polytropic 

model; b) this model. The parameter values are the same of Fig. 4 


obtain X-ray data of the surface brightness as a function of the 
photon energy and position on the sky. Such data could allow 
to derive the gas density, temperature and the abundances of 
several elements as a function of position in the cluster. For the 
energies involved in the X-ray cluster emission the required in- 
strument of these observations must be sensitive to X-ray photon 
energies of at least 7 keV. In relatively nearby clusters, moderate 
spatial and spectral resolutions are necessary and, obviously, the 
sufficient sensitivity to detect the cluster. 


4. Discussion 


It is well known that the hot intracluster gas contains a signifi- 
cant portion of processed gas by stars and shed from galaxies via 
protogalactic winds or via centrally concentrated ram pressure 
stripping. If the second mechanism is predominant, the overall 
abundances of heavy elements could be much less than those 
observed in cluster cores, diminishing the degree of stellar pro- 
cessing required. Protogalactic winds should instead distribute 
metals fairly homogeneously. In this case much processed gas 
is required, with protogalaxies impoverished of perhaps half of 
their initial luminous mass by metal-rich winds. This problem 
regarding the abundance distribution should be clarified by the 
Japanese satellite ASCA (Astro-D) able to obtain spatially re- 
solved X-ray spectral data. Meanwhile, White (1991) analysed 
the possibility to resolve this question evaluating the energet- 
ics associated with metal enrichment in each of these two in- 
jection scenarios. He shows that the consistency of the wind 
scenario depends on the quantity g = €wina/o2 = wfwo,/o2 
where w is the mass fraction of the intracluster gas injected by 
protogalactic winds; fy = 1/3uz,/0% is a factor wich parame- 
terizes the wind velocity,v,,, in terms of a typical stellar veloc- 
ity dispersion o. The most efficient wind-blowing galaxies are 
those relatively early-type galaxies with spheroidal components 
of great luminosity. If gq < 1 the energy deposition by galac- 
tic winds is negligible with respect to cluster collapse. Taking 
w ~ 0.5 and f, < 3 as reasonable rough estimates, the con- 


clusion is that the wind scenario is more plausible in low-o, 
clusters (0, < 500km/s) where a, is comparable to the stellar 
velocity dispersions (0, ~ 300km/s). In high-c, clusters, the 
energy injection by winds is probably negligible, of the order of 
10% also with w = 1. Renzini et al. (1993) claim that it is un- 
likely that a large fraction of the iron mass present in the ICM 
came from ram pressure stripping, essentially because richer 
clusters should show a larger iron mass than poorer clusters due 
to their larger velocity dispersion. An appreciable trend of the 
iron mass-to-light ratio with the cluster richness is not observed. 
Besides, it is necessary to emphasize that the presence of possi- 
ble abundance gradients could be weakened by phenomena of 
subcluster merging. 

As the analysis of the X-ray brightness distributions data 
have shown, the intracluster gas temperature appears to be 
greater than the galaxy temperature (i.e. G < 1). This excess 
of energy seems to be greater in cool clusters ((G si) ~ 0.41) 
than in hot clusters ((Gsi¢) ~ 0.64). The protogalactic winds 
can provide this extra energy to the gas, at difference of the ram- 
pressure stripping mechanism. While this excess of energy can 
be explained in cool clusters or equivalently low velocity dis- 
persions clusters where supernovae-driven protogalactic winds 
can provide this additional energy to the intracluster gas, in hot 
clusters or high velocity dispersion clusters the deposited energy 
by winds is probably negligible as shown by White (1991). 

Some authors (see for example Sarazin 1988; White 1991) 
retain that the average value of 67;¢(~ 0.64), derived by the 
sample of clusters analysed by Jones & Forman (1984), is in 
effect near the unity if one considers that in the isothermal (- 
model, (see Eq. (5)), compares the King approximation instead 
of the isothermal sphere in the description of the cluster poten- 
tial. The King approximation breaks down at large radii where 
varies aS Po, ~ T~> while an isothermal sphere density varies 
aS Ptot ~ 12. This involves that Brit ~ 3/2(Bzit) = 1. If 
so, the hot clusters (high-o,) should not present any energy 
excess of the intracluster medium. Recently, Bahcall & Lu- 
bin (1994) have proposed a similar argument to explain the 
discrepancy between Gspe- and (yj. They assume the actual 
observed average galaxy density distribution in rich clusters : 
par r—2-440.2 for 0. 5h-! <r < 1.5h—'Mpe (Seldner & Pee- 
bles 1977; Peebles 1980) instead of pg ~ r~>. Consequently, 
for r > R, the gas density distribution assumes the following 
slope n ~ pe, mw po 2-4t0.2067i4 — 7—38rit. The correct BF ie 
parameter is related to the parameter (;, determined using the 
King approximation by BF it ~ (1.25 +0.1)Gfi¢ ~ 0.84 40.07 
for Brit ~ 0.67. This value of BF it is compatible with the value 
Bspec ~ 0.94 + 0.08 derived from Lubin & Bahcall (1993). 
Their conclusion is that no significant -discrepancy exists in 
clusters of galaxies and the isothermal G-model provides a con- 
sistent fit to both the X-ray and the optical data. 

However, for r < 0.5h~'Mpc the King model gives an ac- 
curate description of the density profiles. Thus, to provide a 
good fit, the Seldner & Peebles distribution must be normalized 
to give the same density of the King model at r=0.5h~'Mpce. 
When the appropriate normalization factors are considered, a 
smaller correction obtains, and 64;, < 0.8. Besides, it remains 
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difficult to explain with this argumentation the difference be- 
tween groups and rich clusters and in particular to make galaxy 
groups compatible with a @ value ~ 1. The recent analysis of 
FFH on the Coma cluster with the CFF non-polytropic model 
has shown that also without the use of the King approximation 
for the total mass distribution, the value of Gf; is < 1, namely 
Brit = 0.644%. pratically coincident with the value found 
by Hughes et al. (1988), fitting the same X-ray data with the 
isothermal J-model (4; = 0.63 + 0.03). The question of why 
the gas is hotter than the galaxies seems to remain. 

David et al. (1991) examined the dynamics and X-ray evolu- 
tion of the interstellar medium in elliptical galaxies. They found 
that all the models predict Type II supernovae-driven winds at 
early epochs and galactic winds from low-luminosity ellipti- 
cals for most of their lifetime. The evolution of the interstel- 
lar medium in elliptical galaxies has significant effects on the 
observed properties of the intracluster medium. The observed 
iron abundance, oxygen to iron mass ratio and gas tempera- 
ture can all be simultaneously explained if the first generation 
of stars formed with a flat initial mass function. The calcula- 
tions of David et al. (1991) show a rapid decline of the heating 
due to galactic winds from poor to rich clusters. These models 
seem indicate that Type II supernova-driven winds can heat the 
intracluster gas at temperatures above those produced by grav- 
itational heating alone and that consequently the ram-pressure 
stripping is not necessary at significant levels. The questions of 
why 6 < 1 and why there is an increase of its value from poor to 
rich clusters can find an explanation in the simulations of David 
et al. (1991) on the X-ray evolution of the interstellar medium 
in elliptical galaxies with a first generation of stars formed with 
a relatively flat initial mass function. 


5. Conclusions 


The best fitting values of 6 ¢;, obtained from the X-ray brightness 
distributions data, using the isothermal-@ model, range from 
0.41 for groups to 0.64 for rich clusters. From the EXOSAT 
observations, Edge (1989) reports the temperatures for a sample 
of 35 clusters. Two-thirds of the clusters examined by Edge have 
also measured velocity dispersions that give an average value of 
Bspec ~ 0.83. Both the Gr i¢ and Bzpe- values indicate an energy 
excess of the intracluster gas with respect to the energy acquired 
from the gravitational collapse alone. Type II Supernova-driven 
protogalactic winds by early ellipticals with a first generation 
of stars of high masses appear to be the only mechanism able 
to provide this extra energy to the ICM and satisfy the observed 
trend of the ( value increasing from poor to rich clusters (David 
et al. 1991). White (1991) claims that in rich clusters the energy 
contribution to the ICM may be not greater of ~ 10% due to the 
high velocity dispersion with respect to a typical stellar velocity 
dispersion of ~ 300kms~'. 

In this paper we have formulated a static model that includes 
the contribution of galactic winds to the ICM in terms of injec- 
tion of material and heating. This model is a generalization of 
the CFF non-polytropic model used by FFH to jointly fit all 
available X-ray and optical imaging and spectroscopic data of 


the Coma cluster. For the formulation of this model we suppose 
that the binding mass of the cluster is distributed like the galax- 
ies and that the anisotropy of the galaxy velocity distribution is 
constant with the radius. These assumptions seem to be rather 
reasonable. However, to fully explore the consequencies of this 
model, it could be interesting to evaluate the effects of releasing 
both these hypothesis. 

The use of this model for the analysis of spatially and spec- 
trally resolved X-ray surface brightness distributions, as those 
that should be obtained by ASCA (Astro-D), should allow to 
evaluate the consistency of the wind scenario. In particular, the 
model would permit to quantify the contribution of protogalac- 
tic winds to the thermal content of the hot medium in clusters 
of galaxies and verify the importance of the early winds with 
respect to the ram pressure stripping for the metal enrichment 
of the gas. Besides the study of the wind scenario can also con- 
tribute to define the evolution of the member galaxies of the 
clusters. 
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Abstract. We present numerical results of a computer code 
which solves the transfer equation for the synchrotron radiation 
in extragalactic relativistic jets with bends and shock waves. 
We consider for the first time the influence on the emission of 
bent jets by two relativistic effects: aberration of light and time 
delays. We show that the importance of aberration effects de- 
pends on the degree of randomly oriented magnetic field within 
the jet, being null for jets with a completely randomly oriented 
magnetic field. We study the consequences of including time 
delays in bent shocked jets with diverging flow through a three 
dimensional numerical simulation. The time delay produces ro- 
tations and deformations of the shock fronts when viewed in the 
observer’s frame, which result in variations of the length of the 
shocked region and of the apparent relative velocity between 
different fluid lines. Hence the emission, in particular the lin- 
early polarized flux density, associated with the shocked region 
is strongly altered. We present numerical results which com- 
pare those obtained before and after taking into account both 
effects for a particular shocked, helical model jet, and show that 
the corrections are especially important in the polarized flux 
density. 


Key words: galaxies: jets — shock waves — polarization — meth- 
ods: numerical 





1. Introduction 


In Gémez et al. (1993, hereafter Paper I) we described a com- 
puter code which calculates the synchrotron emission from bent, 
shocked relativistic jets. Such jets have been proposed to ex- 
plain phenomena observed on parsec scales in quasars and sim- 
ilar objects. This code considers adiabatic expansion in the jet, 
calculating the evolution of the parameters that determine the 
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synchrotron emission along the jet. Through the assumption of 
a predominantly turbulent magnetic field, the code computes 
the absorption and emission coefficients for the synchrotron 
radiation in cells within the jet, then solves the transfer equa- 
tions along integration columns parallel to the line of sight. 
The code calculates the Stokes parameters of the synchrotron 
radiation along the jet, which, upon convolution with an ap- 
propiate beam (depending on the observing frequency), allow 
us to obtain simulated total and polarized intensity maps. In our 
numerical model we allow the possibility of any form of three 
dimensional geometry for the jet, as well as that of plane perpen- 
dicular shock waves moving along the jet. This provides us with 
the tools to study the influence of bends and shock waves in the 
synchrotron emission from relativistic jets. In Paper I we showed 
that stationary components observed in many extragalactic ra- 
dio sources can be understood as bends in the geometry of the jet 
toward the observer, which produce enhanced emission through 
Doppler boosting. In Gémez et al. (1994, hereafter Paper II) we 
applied our numerical code to the cases of rectilinear and helical 
shocked jets. We presented a multi-epoch and multi-wavelength 
analysis of the application of our numerical code to these hypo- 
thetical sources which can help to understand real sources. We 
then showed that superluminal components can be associated 
with the presence of shock waves, which compress the gas and 
produce enhanced emission. Apparent superluminal motion is 
observed when these relativistic shock waves travel along jets. 


The relativistic effect known as light aberration, responsi- 
ble for the change in the direction of the line of sight between 
different reference systems moving at relativistic relative veloc- 
ities, have been addressed in previous models (Jones et al. 1985; 
Hughes et al. 1989) for the case of rectilinear jets. On the other 
hand, Marscher et al. (1992) showed that another relativistic ef- 
fect present in non-stationary relativistic jets, time delays, can 
be of significant importance in the study of the evolution of the 
emission from shock waves. In an attempt to improve further 
our theoretical-numerical code, we have introduced these two 
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relativistic effects for the case of bent jets, and studied their 
influence on the emission of extragalactic relativistic jets. 


2. Aberration of light 


The synchrotron radiation from relativistic jets depends strongly 
on the orientation at which we observe the jet flow. This orienta- 
tion determines the most important characteristics of the emis- 
sion, such as the total and polarized flux densities, the spectral 
index, etc. Since the jet convects a plasma moving at relativistic 
speeds, we have to consider the Lorentz transformation of the 
direction represented by the line of sight between the observer’s 
frame and the fluid frame, known as aberration of light (see, e.g. 
Rybicki & Lightman 1979). 

Consider a jet fluid element in comoving frame ©’ which 
is moving at a velocity @ (in light units) measured in the ob- 
server’s frame ©. For simplicity, but without loss of generality, 
we consider the velocity ( to lie along the z-axis. If we char- 
acterize the orientation of the line of sight in the © frame by 
the azimuthal angle y, measured from the x-axis and defined as 
positive toward the y-axis, and the elevation angle #, measured 
from the z=0 plane and defined positive toward the positive z- 
axis, the Lorentz transformations between both frames provide 
the orientation of the line of sight as viewed in the fluid frame 
D’ through the following aberration formulae: 


sin @ 





sin 6! = 
a (i — Bcos 6 cos y) ) 





4! V(cos @cos y — 3)?T? + cos? @ sin” y 
cos 0’ = 
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where the primed variables are measured in the fluid frame L’, 
and I is the Lorentz factor of the fluid element in the frame ©. 

We have plotted in Fig. 1, for the case where the azimuthal 
angle yp is zero, the relation between the viewing angles in the 
observer’s and fluid frames for different values of [’. Note that 
there is no loss of generality in Fig. | because aberration only 
affects the angle formed between the line of sight and the direc- 
tion of the velocity vector, which is represented by the angle 0. 
Since the aberration of light is a relativistic effect, the change 
in the direction of the line of sight measured by an observer 
in the frames © and ¥’ depends strongly on the Lorentz factor 
I’. Higher values of [ lead to larger differences between both 
frames. 

Let us consider the consequences of aberration in the syn- 
chrotron emission of relativistic jets. To do this, we estimate 
the magnitude of the corrections that this effect can introduce 
to the results of our numerical model. When solving the trans- 
fer equations for the synchrotron radiation using our numerical 
code, there are two steps where the line of sight appears as a 








40 
6 (°) 
Fig. 1. Relation between the viewing angles in the observer’s (@), and 
fluid (6’) frames, due to aberration effects. The curves correspond to: 


['=10, continuous line; [=8, dashed line; [=6, dashed-dotted line; and 
['=4, dotted line 








parameter: i) in the emission and absorption coefficients, which 
are proportional to B sin v, where B is the magnetic field and 
J is the angle between the line of sight and the magnetic field, 
both measured in the fluid frame (see Paper I); and ii) in the 
computation of the polarization properties, through the angle 
xB, defined as the angle between the direction of the projected 
magnetic field on the plane of the sky, and the axis with respect 
to which we are measuring the polarization angle (see Paper I 
for details). 

The orientation of the magnetic field determines the angles 
% and yg. Therefore, the degree of randomness of the field 
influences the corrections produced by aberration through the 
angles ¥ and yz. If the magnetic field has a completely random 
orientation within the jet, the angles 3 and yg will remain ran- 
dom in orientation and hence, aberration does not introduce any 
correction to the emission in the jet. If the degree of random 
orientation of the magnetic field is substantially reduced, the 
aberration can change drastically the emission properties of the 
jet, becoming a very important effect in the case of an uniform 
magnetic field. 

In order to have a quantitative measure of the corrections in 
the total and polarized flux densities required by the inclusion 
of aberration effects, we have run our numerical code for a rel- 
ativistic jet whose parameters are given in Table | (see Paper 
I for a complete description of every parameter) for two cases: 
i) including aberration, and ii) not including this effect. In both 
cases we have restricted ourselves to the calculation of the emis- 
sion in the plane formed by the line of sight and the ridge line 
of the jet. Since the corrections depend on the degree to which 
the field is ordered, we have performed several sets of paired 
runs for fractions of random field ¢ ranging from 0 to 1.0. Then, 
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Table 1. Physical parameters of the jet models used to quantify the 
corrections by the aberration of light effect 





Value 
0.001 


Parameter 
Noo 


Description 

Relativistic electron density per 
unit energy® (cm~* erg?~') 
Magnetic field strength® (Gauss) 
Cross-sectional radius 

of the jet®* (pc) 

Frequency of observation (GHz) 
Electron energy spectral index 
Velocity of the quiescent flow? 
Opening half-angle of the jet (°) 
Range for the fraction of randomly 
oriented magnetic field 
Exponent in the dependence of 
the magnetic field on radius r 
Redshift of the source 
Rectilinear jet of length 30 pc 
Size of the cells (pc) 

Azimuthal angle (°) 

Elevation angle (°) 





Bo 0.3 
0.1 


22.2 
2.0 
0.98601 
1.0 
0-1 


1.0 


Z 0.6 
Geometry: 








Scells 0.01 
0) 80.0 
h 0.0 





“ At the fiducial point d, and in the fluid frame 
> In the source’s frame and in light units 


for every value of the fraction of random field, we have calcu- 
lated the mean value, for the different positions along the jet, of 
the percentage correction introduced by the aberration of light 
effect in the total and polarized flux densities. 


We have considered a rectilinear jet with a Lorentz factor 
of 6, and a viewing angle in the observer’s frame of 10° for the 
fluid elements in the ridge line of the jet, which corresponds to a 
value of about 90° for the viewing angle in the fluid frame (See 
Fig. 1). That is, when we observe this jet with an angle of 10° 
we are observing the emission that was emitted in a direction 
perpendicular to the velocity of the fluid in the fluid frame. This 
suggests that these corrections can be very important, but as we 
will see later, this depends on the degree of randomness of the 
magnetic field within the jet. 


Figure 2 shows the correction to the total flux density of the 
jet due to relativistic aberration. As can be seen, the corrections 
for high percentages of randomly oriented magnetic fields are 
very small. The corrections are smaller than 2% for ¢ between 
0.9 and 1.0. For a completely random magnetic field (¢ = 1) 
our code estimates a correction of 1.5%, whereas theoretically 
it should be zero. The explanation can be found by taking into 
account that this is the scale of variation of the total flux density 
due to the random behaviour of the magnetic field. In other 
words, this is the intrinsic run-to-run variation of the total flux 
density for a completely randomly oriented magnetic field (for 
the given value of the size of the cells — see Table 1). 

In our numerical applications we usually consider ¢ > 0.7, 
in order to obtain a degree of polarization comparable to the ob- 
servations. For these cases, the corrections due to the aberration 





Total flux density correction (%) 
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Fig. 2. Corrections to the total flux density introduced by the effects of 
relativistic aberration. The models correspond to the parameters listed 
in Table | 


of light do not exceed 10%. However, for smaller percentages 
of randomly oriented field the aberration correction increases 
steeply, reaching values of about 50% for ¢ ~ 0.35. 


We have performed similar calculations for the polarized 
flux density. We show the results in Fig. 3. The polarized flux 
density is extremely sensitive to magnetic field orientation, and 
for a completely randomly oriented magnetic field (¢ = 1) its 
intrinsic indetermination is higher than that for the total flux 
density. Furthermore, in order to be able to compare these cal- 
culations to those, without aberration taken into consideration, 
we have smoothed the data. This procedure introduces a cer- 
tain additional error in the estimation of the corrections. We 
have obtained for the models with ¢ = | an intrinsic fluctua- 
tion (see above) of 6.2%, which includes an error introduced by 
the smoothing of about 3 or 4%. For smaller values of ¢, the 
estimated aberration corrections to the polarized flux increase 
sharply, being about 50% for ¢ = 0.75. Although this correc- 
tion is very important, it is worth noting that the profile of the 
polarized emission along the jet does not change significantly 
as aberration becomes important, but the polarized intensity in- 
creases substantially. For ¢ < 0.4 we find the opposite result, 
that is, the polarized fluxes are comparable in both models but 
the intensity profiles are shifted along the jet axis. If we do not 
consider the effects of aberration, the maximum in the polarized 
flux curve along the jet is located closer to the origin of the jet. 
This behavior can be explained by considering that when the 
randomly oriented magnetic field fraction is small, the factor 
Bsin#v in the absorption and emission coefficients takes on a 
value of ~ B sin 90° if aberration is included, and ~ B sin 10° 
if it is not. Therefore, in the second case the situation is similar 
to the case of a smaller magnetic field (by a factor of sin 10°), 
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Polarized flux density correction (%) 
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Fig. 3. Corrections to the polarized flux density introduced by the ef- 
fects of relativistic aberration. The models correspond to the parameters 
listed in Table | 


which results in a less opaque jet, and the maximum intensity 
position shifts to a position closer to the core of the jet. 

The estimates of the corrections due to aberration we have 
presented thus far should be considered as approximate upper 
limits, because here we have considered a jet model where this 
effect is maximized. 

It should be noted that the percentage of randomly oriented 
magnetic field can be significantly smaller in regions excited 
by shock waves than in the rest of the jet, and therefore greater 
corrections due to aberration in these regions are to be expected. 
Nevertheless, the presence of a shock wave does not necessar- 
ily imply a reordering of the magnetic field: this depends on 
the initial percentage of randomly oriented magnetic field and 
on the strength of the shock wave. For instance, if a jet has a 
relatively small fraction of randomly oriented magnetic field, 
say 30-40%, and a uniform component directed at right angles 
to the shock front, a shock wave will introduce a higher degree 
of randomization, instead of reducing it. 


3. Time delays 


The time delay effect is relevant when we consider relativistic 
jets which are not stationary with time, for example when a 
shock wave, or any other kind of inhomogeneity, travels along 
the jet. In such cases, the difference in time of arrival at the 
observer of emission from different positions in the jet will 
produce changes in the shape and width of shocked regions, 
as viewed in the observer’s frame, as a function of the viewing 
angle. 

Let us consider firstly a rectilinear jet viewed at a certain 
angle @, defined as the angle made between the line of sight and 
the ridge line of the jet. Let us assume that a perturbation, e.g. a 


shock wave, is travelling along the ridge line of the jet between 
positions | and 2 at a speed (c. The difference in arrival times in 
the observer’s frame, AT’, of the radiation emitted at positions 
1 and 2, is equal to that measured in the source’s frame, At, 
reduced by the time required by the radiation to propagate the 
distance GcAt cos 6. Additionally, we need to include the time 
dilation between both frames due to the source redshift, which 
is a multiplicative factor of (1+z). Thus we have 
AT = At(1 — Bcos6)(1 + z) (5) 

Because of this difference in time between two events in the 
observer’s and source’s frame, the shape of a feature, such as a 
shocked area, may appear different in the two frames, and is a 
function of the viewing angle 0. 

We could interpret the difference in times through a differ- 
ence in velocities between the two reference systems. Since the 
perturbation travels along a distance JcAt in a period of time 
AT in the observer’s frame, the observer would measure that 
the perturbation moves with a velocity 


_ Bcht | Bec 
~ AT” (1—Bcosé\(1 +z) 


V 





(6) 


which differs by a factor (I—@ cos 0)~'(1+z)~! from the velocity 
measured in the source’s frame, v = (Gc. Since the jet convects 
divergent flow, different fluid lines have different viewing an- 
gles, and therefore different velocities in the observer’s frame. 
For instance, if we consider a shock wave travelling along the 
jet, the time delay will produce a rotation of the two (forward 
and reverse) shock fronts as viewed in the observer’s frame, 
as the shock advances along a relativistic jet (Marscher et al. 
1992). 

We will center our attention on the case of a shock wave 
travelling along a bent jet, without any restriction to the three 
dimensional geometry of such a jet. In order to apply Eq. (6) to 
bent jets we consider that the orientation of the jet with respect to 
the observer depends on the position along the jet for every fluid 
line. Hence the viewing angle 6; for a fluid element travelling 
along a fixed fluid line in the jet is a function of time 6¢;(T). 
Then we can apply Eq. (6) to determine the shape of the shocked 
region as a function of time in the observer’s reference frame. 
Since the shocked region is determined by the portion of gas 
situated between the forward and the reverse shock fronts, we 
only need to know the position of both shock fronts as a function 
of time. As most shocks in diverging flows are self-similar, the 
velocity of the reverse shock in the source’s frame will be 2G.) 
where (,;, is the velocity of the forward shock front, and x a 
certain constant. Assuming that the shocked region has an initial 
length /;;,. in the source’s frame, the position of both fronts 
along the jet in the observer’s reference frame can be calculated 
as a function of the observer’s time. Integrating along each fluid 
line, we obtain 


pf (T) i [ Bshe te 
ft 9 (1 — Bsn cos Op,(7)\(1 + 2) 
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LBsnc dr 
(1 —2BsncosOs(T))1+z) 
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where 


(8) 


lsho(1 — Bsn cos @¢,(0))(1 + z) 
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Here D},(T) and D’;,(T’) are the positions along the jet of the 
forward and reverse fronts, respectively, at a given time T’ and 
for a given fluid line. The position of the reverse front is obtained 
by integrating up to a certain time 7}. which is the time of obser- 
vation 7’ reduced by the time — measured in the observer’s frame 
— needed by the forward shock to travel a distance equal to the 
initial length of the shocked region, following a rectilinear path 
determined by the initial orientation @;(0) for each fluid line. 
Note that we have assumed that the initial length /,,,, and G,;, 
are the same in the source’s frame for every fluid line. With this 
definition for T;., the initial length of the shocked region for ev- 
ery fluid line in the observer's frame can be different from |p, 
since the time delay affects the forward shock during the time 
necessary for the reverse shock front to cross the fiducial point 
of the jet, 7’ — T;, (note that this depends also on the particular 
fluid line in the jet). 

Since jets convect divergent flow, time delays depend on the 
fluid line, and therefore a definition of apparent velocity of the 
shocked region as a whole is not possible at all. Nevertheless, 
considering that the VLBI observations do not have resolution 
sufficient to observe these effects, we could define a mean ap- 
parent velocity of the component associated with the shocked 
region. 

Hence time delay corrections in bent jets require the knowl- 
edge of the angle between the instantaneous flow velocity vector 
and the line of sight @ ¢,(T’) at every time and for every fluid line 
in the jet, which allows the solution of Eqs. (7) and (8). 


3.1. A three dimensional numerical solution 


We can define the three dimensional geometry of any jet through 
a set of points in space which define the ridge line of the 
jet. Each of these points are characterized by their coordinates 
(Zri, Yrl, Zr) in the source’s reference system. The number of 
points required to describe the geometry depends on the geome- 
try itself, that is, for rectilinear jets we need only two points, but 
if the jet bends strongly we need more points to define accurately 
the curvature. 

Consider in such a reference frame a unit vector s;, which 
represents the direction of the line of sight (see Fig. 4). This 
vector can be specified by two angles: ¢;,, measured in the 
plane z=0 from the x-axis and positive toward the y-axis; and 
W1s, the angle between s;, and the plane z=0, positive toward the 
z-axis. In order to account for time delay through Eqs. (7) and 
(8), we need to calculate numerically the angle 6 r;(T), defined 
as the angle between the vector s;, and the velocity vector of 
each fluid element. 

Let us consider first the case of the ridge line. We define the 
orientation of the velocity vector for each point in the ridge line 


Zz 4 
z'4 














Fig. 4. Reference systems and definition of fluid lines for the three 
dimensional numerical solution 


associated with a fluid element, v? , aS a unit vector that con- 
nects the point (7) with the next point (7 + 1), calculated through 
the coordinates of the two points. Knowing the orientation of 
the vectors v? and s;,, we can calculate the angle G for each 
point, and apply Eqs. (7) and (8) to obtain the positions along 
the ridge line of both the forward and the reverse fronts of the 
shock wave, at a given time 7’. Note that the index (2) can be 
associated with a certain time T. 

Let us consider now fluid lines which are not on the ridge 
line of the jet. Consider as a reference a point (7) in the ridge 
line of the jet. We can define a new reference system at this 
point through a rotation of an angle ¢? along the z-axis of the 
source’s system and a rotation of an angle y along the y’-axis, 
resulting from the first rotation along the z-axis of the source’s 
system (see Fig. 4). In this new reference system the fluid lines 
can be determined by an angle ar; measured in the plane z’ = 0, 
and an angle y+, with respect to the x’-axis. To account for all 
the possible fluid lines in the jet, the angle af; will vary from 0 
to 27, and y , will vary from 0 to ¢;, the half-opening angle of 
the jet. 

To solve Eqs. (7) and (8) for a given fluid line characterized 
by the values af; and yy), we define a unit vector in the new 
reference system by the coordinates 


x’ =cos yi (10) 


(11) 


y’ = sin pf COs af; 


2’ =sinyys sinas (12) 





